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INTRODUCTION:  The  subject  of  this  research  is  the  need  for  improved  osteoporosis  therapies.  The 
purpose  of  this  study  is  to  test  the  efficacy  of  parathyroid  hormone  (PTH)  mono-therapy  and  PTH  +  anti- 
cataholic  combination  therapies  on  Nmp4~/~  and  wild  type  (WT)  mice.  We  have  previously  determined  that 
Nmp4  represses  the  response  of  hone  to  osteoanaholics^"^.  The  scope  of  this  research  comprises  the 
following  specific  aims:  (i  &  ii)  to  determine  the  impact  of  Nmp4  on  the  efficacy  of  PTH  mono-  and 
comhination  therapies  with  hisphosphonates  and  the  selective  estrogen  receptor  modulator  (SERM)  raloxifene 
in  ovariectomized  (ovx)  mice;  (iii)  to  determine  the  cell  type-specific  contributions  to  the  enhanced  response  of 
the  Nmp4~/~  mouse  to  these  osteoporosis  therapies. 

BODY: 

During  Year  2  we  successfully  accomplished  a  number  of  objectives  in  the  following  Tasks  and  Subtasks: 

i.  Task  2:  Conduct  PTH  mono-therapy  and  PTH  combination  therapies  with  hisphosphonates:  A 
revised  manuscript  describing  the  impact  of  Nmp4  on  ovx-induced  bone  loss  and  the  response  of 
ovx  mice  to  PTH  mono-therapy  is  in  review  at  Molecular  Endocrinology  (see  Appendix  Manuscript 
Number  ME-14-1406R1).  Note:  this  manuscript  contains  all  details  for  the  materials  and  methods 
used  to  generate  the  PTH+anti-catabolic  data  described  below. 

ii.  Task  3:  Conduct  PTH  combination  therapy  with  raloxifene:  The  treatment  of  the  WT  and  NmpyT 
mice  with  the  PTH  +  anti-catabolic  therapies  (raloxifene  and  hisphosphonates)  is  finished  and  the 
analyses  underway.  The  PTH+raloxifene  combination  therapy  has  been  identified  as  the  optimum 
combination  treatment  to  be  used  in  Task  4  (see  below  for  details  of  results  and  ongoing  analyses). 

iii.  Subtasks  4.2. a  and  4.3.a:  Breed  54  female  Nmp4fl/fl3.6CoTCre+  &  54  Nmp4fl/fl-Cre-  mice:  We 
are  backcrossing  our  Nmp4®°Yflox  mice  onto  a  C57BL/6J  for  production  of  conditional  knockout 
animals  (see  below  for  details). 

iv.  Subtask  6.2:  Culture,  expand,  and  characterize  the  phenotype  of  both  WT  and  NmpyC 
mesenchymal  stem/progenitor  cells  (MSPCs):  These  experiments  are  described  in  the  revised 
manuscript  ME-14-1406R1  (Appendix)  and  we  have  started  experiments  designed  to  characterize 
the  impact  of  raloxifene  on  MSPCs  (see  below) 

Overview  of  Objeetives  and  Results  During  Year  2: 

During  Year  1  we  determined  that  the  exaggerated  response  to  anabolic  doses  of  PTH  was  preserved  in  ovx 
Nmp4~/~  mice  but  that  disabling  this  transcription  factor  did  not  protect  them  from  ovx-induced  bone  loss 
without  therapy.  During  Year  2  we  addressed  whether  various  PTH+anti-catabolic  combination  therapies  are 
more  efficacious  in  osteoporotic  than  wild  type  (WT)  mice.  Animals  were  ovx  at  i2wks.  Therapies  were 

initiated  at  i6wks  of  age  including  PTH  only,  alendronate  (ALN)  only,  raloxifene  (RAL)  only,  zolendronate 
(ZOL)  only,  and  the  combination  treatments  PTH + ALN,  PTH + RAL,  PTH+ZOL.  The  vehicle  control  (VEH)  was 
comprised  of  all  the  carriers  used  in  the  study.  At  24wks  of  age  (8wks  treatment)  mice  were  euthanized  for 
analysis. 

Doses: 

•  PTH:  3opg/kg/d 

•  ALN:  ipg/kg/d 

•  ZOL:  8opg/kg  IX  dose: 

_ RAL:  img/kg/d 


Doses  of  anti-resorptive  agents  are  based  on  human  clinical  doses. 

The  standard  ALN  dose  for  treatment  of  osteoporosis  typically  given  as  either  a  daily  (lomg)  or  weekly  (yomg) 
dose.  Based  on  a  60kg  individual  this  is  roughly  1.17  mg/kg/ week.  The  human  dose  is  oral  and  has  an 
estimated  bioavailability  of  around  0.6%,  meaning  that  the  absorbed  dose  is  roughly  0.007  mg/kg/week  (or  7 
pg/kg/week).  We  dosed  via  injection,  assuming  100%  absorption,  thus  we  delivered  ALN  at  ipg/kg/day^-^. 

RAL  is  typically  given  as  a  60  mg  daily  dose.  Based  on  a  60  kg  patient,  the  dose  would  be  1  mg/kg/day.  The 
assumption  is  100%  absorption  thus  the  full  dose  is  used  when  injectings-^. 

ZOL  is  typically  given  yearly  at  a  dose  of  5  mg.  Based  on  a  60  kg  patient,  the  dose  is  0.083  mg/kg.  Our  single 
dose  of  80  Ml/%  approximates  this  amount^’  9. _ 
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SUMMARY  OF  KEY  FINDINGS  FOR  YEAR  2: 

We  have  generated  an  extensive  data  set  from  the  combination  experiments  and  although  the  analysis  is  not 
complete  the  preliminary  statistical  evaluation  provides  an  exciting  emerging  picture,  which  is  presented  here. 
Note  that  Nmp4  targets  cancellous  bone  thus  we  will  specify  what  skeletal  compartments  comprise  the 
endpoint  under  consideration. 

•  The  PTH+RAL  combination  therapy  typically  added  more  bone  to  both  WT  and  Nmp4~/~  ovx  skeletons 
than  the  other  treatments  (whole  body  [WB]  and  spine  BMD  and  femoral  and  spine  BV /TV) 

•  The  PTH+RAL  treatment  was  the  most  effective  at  out-performing  the  PTH  mono-therapy  in  both 
genotypes  (whole  body  [WB]  and  spine  BMD  and  femoral  and  spine  BV/TV) 

•  Disabling  Nmp4  significantly  improved  the  efficacy  of  the  PTH+RAL  therapy  but  typically  had  little  to 
no  effect  on  the  PTH+bisphosphonate  treatments  (whole  body  [WB]  and  spine  BMD  and  femoral  and 
spine  BV/TV) 

•  The  anti-catabolic  mono-therapies,  as  expected,  added  only  a  modest  amount  of  bone  to  the  ovx  WT 
and  Nmp4'/~  skeletons,  compared  to  PTH  mono-therapy.  Disabling  Nmp4  did  not  enhance  the  response 
of  the  ovx  skeleton  to  the  anti-catabolic  mono-therapies. 

•  The  Nmp4~/~  PTH+RAL  treatment  group  exhibited  the  largest  pool  of  bone  marrow  osteoprogenitors 
[identified  by  CD457CDi05+/Nestin+/CDi46+]. 

•  Raloxifene  enhanced  the  mineralization  capacity  of  WT  mesenchymal  stem/progenitor  cells  (MSPCs) 
but  did  not  further  augment  the  precocious  Nmp4~/~  MSPC  mineralization. 

STATSTICAL  TREATMENT; 

We  have  tested  7  therapies  and  a  vehicle  control  using  two  genotypes  of  mice  yielding  a  total  of  16  treatment 
groups.  To  date  we  have  used  the  following  statistical  methods: 

>  Raw  data: 

♦  1.1CT  trabecular  data:  Cancellous  bone  is  the  primary  target  of  Nmp4.  To  determine  if 
there  was  a  genotype  x  treatment  interaction  for  any  of  the  combination  treatments  we 
performed  a  series  of  2W  ANOVAs  using  Genotype  and  Treatment  (veh  vs  PTH+anti- 
catabolic;  PTH  vs  PTH+anti-catabolic)  as  the  independent  variables  and  femoral  BV/TV, 
femoral  Tb  N  and  Tb  Th,  or  spine  BV/TV  as  the  dependent  endpoints.  We  set  p<o.oi  to 
avoid  Type  I  errors.  Similarly,  to  determine  if  the  combination  treatments  provided  a 
synergistic  effect  over  their  respective  mono-therapies  we  performed  a  series  of  2W  ANOVAs 
using  PTH  and  the  anti-catabolic  treatments  as  the  independent  variables. 

♦  HOT  cortical  data:  There  was  no  indication  of  a  genotype  x  treatment  interaction  therefore 
the  raw  cortical  data  was  treated  with  a  iW  ANOVA  followed  by  a  Tukey-Kramer  HSD  post 
hoc  test  to  compare  means. 

>  %  Change  BMD  data:  For  these  data  we  used  the  Kruskal-Wallis  iW  ANOVA  by  ranks  for  non- 
parametric  data. 

>  FACS  data:  The  raw  FACS  data  is  reported  as  %  of  total  cells  analyzed.  We  used  the  Kruskal- 
Wallis  iW  ANOVA  followed  by  nonparametric  comparisons  with  control  using  the  Dunn  method  for 
joint  ranking. 
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KEY  FINDING:  Mice  receiving  PTH+RAL  and  PTH+ZOL  showed  the  largest  BMD  increase  after  ovx. 


A.  Percent  Change  Whole  Body  (WB)  BMD 


B.  Percent  Change  Spine  BMD  |  |  wt 


Figure  i:  Percent  change  [A]  whole  body  bone  mineral  density  (WB  BMD)  (i6wks  to  24wks  of  age)  and  [B]  spine  BMD  (L3-L5)  for  WT 
and  Nmp4~/~  [KO]  mice  under  various  treatment  groups.  The  PTH+RAL  and  PTH+ZOL  stimulated  the  largest  changes  in  WB  and  spine 
BMD.  The  data  were  analyzed  using  the  Kruskal- Wallis  test  for  non-parametric  statistics  (p<  o.oooi)  followed  by  a  comparison  for  all 
pairs  using  Dunn  method  for  joint  ranking  post  hoc  test.  The  data  represents  average±SD,  n=7-i2  mice/group. 


A.  Whole  Body  BMD  24wks 


B.  1W  ANOVA 


/Comparison  of 

all  pairs 
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A 
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B  C 
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B  C 

WT-PTH 

C  D 

KO-RAL 

C  D  E 

WT-ZOL 

CDE 

KO-ZOL 

D  E 

KO-ALN 

D  E 

WT-RAL 

D  E 

WT-AIN 

E 

KO-VEH 

E 

WT-VEH 

E 

Levels  not  connected  by  the  same 
letter  are  significantly  different 


C.  Spine  BMD  24wks 


D.  1W  ANOVA 
/Comparison  of  all  pairs 

KO  PTH.RAL  A 
WT-PTHrRAt  A 
KO  PTH«ZOL  A  B 
WT  PTHtZOL  A  B 
KO-PTH  B  C 

KO  PTH+ALN  B  C 
WTPTH»ALN  C  D 
WT-PTH  COE 

KO-RAL  D  E  F 

WT-ZOL  D  E  F  G 

WT-RAL  0  E  F  G 

KO-ZOL  D  E  F  G 

KO-VEH  E  F  G 

KO-ALN  E  F  G 

WT-ALN  F  G 

WT-VEH  G 

Levels  not  connected  by  the  same 
letter  are  significantly  different 


Figure  1  illustrates  the  percent 
change  of  WB  and  spine  BMD  of 
the  WT  and  Nmp4~/~  mice  from  the 
16  different  treatment  groups.  The 
use  of  the  Kruskal- Wallis  test  for 
non-parametric  statistics  showed  a 
significant  difference  between  the 
groups  (p<o.oooi).  The  use  of  a 
Dunn’s  post-hoc  test  comparing  all 
pair-wise  comparisons  indicated 
that  the  PTH+RAL  and  PTH+ZOL 
groups  were  distinct  in  that  they 
were  different  from  most  of  the 
other  treatments  (statistical  output 
not  shown  in  graphs  for  simplicity). 

Figure  2  compares  the  raw  WB  and 
spine  BMD  of  the  different 
treatment  groups  at  the  end  of  the 
experiment  (8wks  of  therapy, 

24wks  of  age).  The  use  of  a  iW 
ANOVA  followed  by  a  Tukey- 
Kramer  HSD  shows  that  the 
PTH+RAL  and  the  PTH+ZOL 
treatments  were  the  most  effective 
at  adding  bone  after  ovx. 


Figure  2:  [A]  WB  BMD  (24wks  of  age)  and  [B]  spine  BMD  (L3-L5)  for  WT  and  Nmp4~/~ 
[KO]  mice  under  various  treatment  groups.  [C],  [D]  The  data  were  analyzed  using  a  iW 
ANOVA  followed  by  a  comparison  for  all  pairs  using  a  Tukey-Kramer  (p<0.05).  The  mice 
from  the  PTH+RAL  (and  PTH+ZOL)  exhibited  the  largest  WB  and  spine  BMDs.  The  WT 
mice  exhibited  a  S5mergistic  response  to  the  PTH+RAL  combination  therapy  (see  Table  1, 
next  page).  The  data  represents  average±SD,  n=7-i2  mice/group. 


These  results  raise  the  question  as 
to  whether  any  of  the  combination 
treatments  show  a  synergistic 
response  at  the  levels  of  WB  and 
spine  BMD. 


6 


KEY  FINDING:  The  PTH+RAL  therapy  was  the  only  treatment  to  exhibit  a  synergistic  effect  on  total  WB 
and  spine  BMD 


TABLE  1 


WHOLE  BODY  BMD 

THERAPY 

p-value  PTH  Treatment 

p-value  Anti-catabolic 
Treatment 

p-value  PTH  x  Anti- 
catabolic  interaction 

PTH+ALN  [WT  mice] 

<0.0001 

0.04 

0.65 

PTH+ALN  [NmpT^^  mice] 

<0.0001 

0.50 

0.64 

PTH+RAL  [WT  mice] 

<0.0001 

<0.0001 

0.0007 

PTH+RAL  [NmpT^^  mice] 

<0.0001 

<0.0001 

0.1025 

PTH+ZOL  [WT  mice] 

<0.0001 

0.0003 

0.8258 

PTH+ZOL  [Nmp4'^'  mice] 

<0.0001 

0.0023 

0.3662 

SPINE  BMD  (L3-L5) 

THERAPY 

p-value  PTH  Treatment 

p-value  Anti-catabolic 
Treatment 

p-value  PTH  x  Anti- 
catabolic  interaction 

PTH+ALN  [WT  mice] 

<0.0001 

0.02 

0.94 

PTH+ALN  [NmpT^^  mice] 

<0.0001 

0.69 

0.79 

PTH+RAL  [WT  mice] 

<0.0001 

<0.0001 

0.003 

PTH+RAL  \Nmp4^^  mice] 

<0.0001 

<0.0001 

0.0175 

PTH+ZOL  [WT  mice] 

<0.0001 

<0.0001 

0.3807 

PTH+ZOL  \Nmp4^'  mice] 

<0.0001 

0.0101 

0.0561 

To  determine  if  the  combination  treatments  provided  a  synergistic  effect  over  their  respective  mono-therapies 
we  performed  a  series  of  2W  ANOVAs  using  PTH  and  the  anti-catabolic  treatments  as  the  independent 
variables  and  either  WB  BMD  or  spine  BMD  (at  8wks  of  therapy,  24wks  of  age)  as  the  dependent  endpoint  (see 
Table  1).  Note  that  these  analyses  used  the  actual  BMD  value  and  not  the  %  change.  We  set  p<o.oi  to  avoid 
Type  I  errors.  Only  the  PTH+RAL  combination  therapy  yielded  a  synergistic  effect  for  the  WT  (WB  BMD/spine 
BMD).  The  Nmp4'/~  showed  a  nearly  significant  synergistic  effect  for  spine  BMD. 

The  PTH+  bisphosphonate  therapies  (ALN  and  ZOL)  did  not  act  synergistically  on  either  the  WT  and  Nmp4~/~ 
mice.  ALN  induced  the  weakest  response  in  both  genotypes  for  these  endpoints  (%  change  BMD  and  BMD, 
Figures,  1,  2,  and  Table  1). 

Potential  interpretations  of  some  of  the  data: 

•  WB  BMD  is  primarily  cortical  bone  (Nmp4  non-target)  whereas  the  spine  BMD  (L3-L5)  has  a 
comparatively  larger  proportion  of  cancellous  bone  (Nmp4  target).  This  may  explain  the  lack  of  a 
synergistic  effect  for  the  Nmp4~/~  mice  for  the  WB  BMD  parameter  but  near  synergistic  effect  for  spine 
BMD. 

•  The  difference  in  the  response  to  ALN  mono-therapy  between  WT  and  Nmp4~/~  mice  is  consistent  with 
our  previous  observation  that  null  osteoclasts  are  more  active  than  WT  cells^.  These  differences  may  be 
overwhelmed  in  the  presence  of  the  stronger  bisphosphonate  ZOL. 
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KEY  FINDING:  Disabling  Nmpg  enhanced  PTH+RAL-induced  increases  in  femoral  BV/TV  and  increased 
the  synergistic  action  of  PTH  with  RAL  and  ZOL  in  both  the  femur  and  spine. 


A.  Femoral  BV/TV  24wks 


B.  L5  BV/TV  24wks 


WT 

Nmp4 


-/- 


PTH-fZQL  PTH+RAL 
'KO  synergy  *GxT 


C.  Femoral  Tb  N  24wks 


PTHfRAL 

*GxT 


PTm-ZOL  PTH4-RAL 
*KO  synergy  'KO  synergy 


PTH*Z0L 
'KO  synergy 


Figure  3:  [A]  Femoral  BV/TV  (24wks  of  age);  [B]  spine  BV/TV  (L5)  [C]  Femoral  trabecular  number  (Tb  N);  and  [D]  Femoral 
trabecular  thickness  (Tb  Th)  for  WT  and  Nmpp/'  mice  under  various  treatment  groups.  The  data  were  analyzed  using  a  series  of  2W 
ANOVAs.  The  PTH  combination  treatment  groups  were  compared  for  genotype  x  treatment  interactions  (G  x  T).  They  were  also 
compared  for  treatment  synergy  (PTH  x  anti-catabolic)  for  the  WT  and  Nmp4~/~  mice  separately.  The  detailed  results  of  the  statistical 
analyses  are  shown  in  Tables  2-4.  There  was  a  G  x  T  interaction  in  the  PTH+RAL  therapy  (femoral  BV/TV  and  Tb  N,  veh  vs  PTH+RAL). 
The  Nmp4'/'  mice  exhibited  a  s}mergistic  response  to  the  PTH+RAL  therapy  (femoral  BV/TV,  femoral  Tb  N,  L5  BV/TV)  and  to  the 
PTH+ZOL  therapy  (femoral  BV/TV,  femoral  Tb  Th,  spine  BV/TV).  The  WT  mice  showed  a  synergistic  response  to  PTH+RAL  in  the 
spine.  The  data  represents  average±SD,  n=7-i2  mice/group. 


Table  21  FEMUR  BV/TV 


Therapy 

Independent  variable  #1 

Independent  variable  #2 

Interaetion  var#i  x  var#2 

PTH+RAL:  Gene  x  TreaV 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=o.ooio 

PTH+RAL:  Gene  x  Treat'’ 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=o. 048,8 

PTH+RAL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<o.oooi 

PTH  X  RAL:  p=o.ooio 

PTH+RAL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<o.oooi 

PTH  X  RAL:  p=o.Oi39 

PTH+ZOL:  Gene  x  TreaV 

Genotype:  p=0. 00,89 

Treatment:  p<o.oooi 

GxT:  p=o.076o 

PTH+ZOL:  Gene  x  Treat'’ 

Genotype:  p=0.0059 

Treatment:  p<o.oooi 

GxT:  p=o.9ioi 

PTH+ZOL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<o.oooi 

PTH  X  ZOL:  p=o.oo67 

PTH+ZOL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<0.000i 

PTH  X  ZOL:  p=o.oi88 

PTH+ALN:  Gene  x  Treat” 

Genotype:  p=o. 029,8 

Treatment:  p<o.oooi 

GxT:  p=o.7,344 

PTH+ALN:  Gene  x  Treat'’ 

Genotype:  p<o.oooi 

Treatment:  p=0.0779 

GxT:  p=o.oi97 

PTH+ALN:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ALN  treatment:  p=0.5S40 

PTH  X  ALN:  p=o.9i48 

PTH+ALN:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ALN  treatment:  p=0.000i 

PTH  X  ALN:  p=o.i944 

a:  Genot)4)e  x  treatment  comparing  vehicle-treated  mice  to  PTH+  anti-catabolic  therapy 
b:  Genot)4)e  x  treatment  comparing  PTH-treated  mice  to  PTH+  anti-catabolic  therapy 
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Table  3;  FEMUR  Tb  N  and  Tb  Th 


Therapy 

Independent  variable  #1 

Independent  variable  #2 

Interaetion  var#i  x  var#2 

Femur  Tb  N 

PTH+RAL:  Gene  x  TreaH 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=0.0090 

PTH+RAL:  Gene  x  Treat^ 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=o.i641 

PTH+RAL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<0.000i 

PTH  X  RAL:  p<o.oooi 

PTH+RAL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<0.000i 

PTH  X  RAL:  p=0.03i3 

PTH+ZOL:  Gene  x  TreaH 

Genotype:  p=0.0046 

Treatment:  p<o.oooi 

GxT:  p=o.404i 

PTH+ZOL:  Gene  x  TreaF 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=o.4838 

PTH+ZOL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p=0.0030 

PTH  X  ZOL:  p=o.i69i 

PTH+ZOL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p=o.ooo4 

PTH  X  ZOL:  p=o.i4i5 

Femur  Tb  Th 

PTH+RAL:  Gene  x  TreaH 

Genotype:  p=0.0337 

Treatment:  p<o.oooi 

GxT:  p=o.0422 

PTH+RAL:  Gene  x  TreaF 

Genotype:  p<0.000i 

Treatment:  p<o.oooi 

GxT:  p=o.7730 

PTH+RAL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<0.000i 

PTH  X  RAL:  p=0.3820 

PTH+RAL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<0.000i 

PTH  X  RAL:  p=0.6i40 

PTH+ZOL:  Gene  x  TreaH 

Genotype:  p=o.oi27 

Treatment:  p<o.oooi 

GxT:  p=o.oi66 

PTH+ZOL:  Gene  x  TreaF 

Genotype:  p<0.000i 

Treatment:  p<o.oooi 

GxT:  p=O..S4.S8 

PTH+ZOL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<o.oooi 

PTH  X  ZOL:  p=o.ooo2 

PTH+ZOL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<o.oooi 

PTH  X  ZOL:  p=o.oi72 

a:  Genot)?pe  x  treatment  comparing  vehicle-treated  mice  to  PTH+  anti-catabolic  therapy 
b:  Genot)qDe  x  treatment  comparing  PTH-treated  mice  to  PTH+  anti-catabolic  therapy 


Table  4:  L5  BV/TV 


Therapy 

Independent  variable  #1 

Independent  variable  #2 

Interaetion  var#i  x  var#2 

PTH+RAL:  Gene  x  TreaH 

Genotype:  p=0.0002 

Treatment:  p<o.oooi 

GxT:  p=o.o.04i7 

PTH+RAL:  Gene  x  Treat^ 

Genotype:  p<o.oooi 

Treatment:  p<o.oooi 

GxT:  p=o.5999 

PTH+RAL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<o.oooi 

PTH  X  RAL:  p<o.oooi 

PTH+RAL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

RAL  treatment:  p<o.oooi 

PTH  X  RAL:  p=o.ooo2 

PTH+ZOL:  Gene  x  TreaH 

Genotype:  p=0.0205 

Treatment:  p<o.oooi 

GxT:  p=o.25,34 

PTH+ZOL:  Gene  x  TreaF 

Genotype:  p=0.0026 

Treatment:  p<o.oooi 

GxT:  p=o.8i7.S 

PTH+ZOL:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<0.000i 

PTH  X  ZOL:  p<o.oooi 

PTH+ZOL:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ZOL  treatment:  p<o.oooi 

PTH  X  ZOL:  p=0.0235 

PTH+ALN:  Gene  x  TreaH 

Genotype:  p=0.0045 

Treatment:  p<o.oooi 

GxT:  p=0.4050 

PTH+ALN:  Gene  x  Treat*" 

Genotype:  p=0.0002 

Treatment:  p=o.i28i 

GxT:  p=o.5i25 

PTH+ALN:  KO  Synergy 

PTH  treatment:  p<o.oooi 

ALN  treatment:  p=o.8438 

PTH  X  ALN:  p=0.3830 

PTH+ALN:  WT  Synergy 

PTH  treatment:  p<o.oooi 

ALN  treatment:  p=0.0472 

PTH  X  ALN:  p=o.3i3i 

a:  Genot)?pe  x  treatment  comparing  vehicle-treated  mice  to  PTH+  anti-catabolic  therapy 
b:  Genot)qDe  x  treatment  comparing  PTH-treated  mice  to  PTH+  anti-catabolic  therapy 


•  To  determine  if  the  combination  treatments  provided  a  synergistic  effect  over  their  respective  mono¬ 
therapies  we  performed  a  series  of  2W  ANOVAs  using  PTH  and  the  anti-catabolic  treatments  as  the 
independent  variables  and  Femoral  BV/TV  (Table  2),  Femoral  Tb  N  and  Tb  Th  (Table  3),  or  Spine  BV/TV 
(Table  4)  as  the  dependent  endpoints  (8wks  of  therapy,  24wks  of  age).  We  set  p<o.oi  to  avoid  Type  I  errors. 

>  Nmp4~/'  mice  show  a  synergistic  response  to  PTH-fRAL  (Femoral  BV/TV,  Femoral  Tb  N,  and  Spine 
BV/TV) 

>  The  WT  mice  show  a  synergistic  response  to  PTH-fRAL  (Spine  BV/TV) 

>  Nmp4~/'  mice  show  a  synergistic  response  to  PTH-fZOL  (Femoral  BV/TV,  Femoral  Tb  Th,  and  Spine 
BV/TV) 

•  To  determine  if  there  was  a  genotype  x  treatment  interaction  for  any  of  the  combination  treatments  we 
performed  a  series  of  2W  ANOVAs  using  Genotype  and  Treatment  (veh  vs  PTH -1- anti-catabolic;  PTH  vs 
PTH -I- anti-catabolic)  as  the  independent  variables  and  Femoral  BV/TV  (Table  2),  Femoral  Tb  N  and  Tb  Th 
(Table  3),  or  Spine  BV/TV  (Table  4)  as  the  dependent  endpoints.  We  set  p<o.oi  to  avoid  Type  I  errors. 

>  Genotype  x  treatment  interactions  were  observed  for  femoral  BV/TV  and  femoral  Tb  N  under  the 
PTH-i-RAL  therapy  when  compared  to  the  vehicle-treated  mice 

>  No  genotype  x  treatment  interactions  were  obtained  for  the  PTH-i-ZOL  therapy 
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KEY  FINDING:  Disabling  Nmpg  enhanced  the  distinction  between  the  PTH+RAL-induced  increase  in 
cortical  bone  area  and  the  other  combination  treatments 


A.  Cortical  Bone  Area  B.  1W  ANOVA 

/Comparison  of  all  pairs 

KO  PTH+RAL  A 
WT  PTW+RAL  A  B 
KO  PTH+ALN  ABC 
KO  PTH+ZOL  ABC 
WT  PTH+ZOL  A  B  C  D 
WT  PTH+ALN  A  B  C  D 
KO  P  TH  B  C  D  E 

WT  P  TH  B  C  D  E 

KO  RAL  C  D  E 

WTALN  D  E 

WT  ZOL  D  E 

WTRAL  0  E 

KO  ZOL  D  E 

WTVEH  E 

KO  ALN  E 

KOVEH  E 

Levels  not  connected  by  the  same 
letter  are  significantly  different 

Figure  4:  [A]  Femoral  cortical  bone  area  of  ovx  WT  and  Nmp4'/~  mice  at  24wks  of  age,  8wks  of  treatment.  Statistical  analysis  included 
a  iW  ANOVA  (p<  o.oooi)  followed  by  a  Tukey-Kramer  HSD  post-hoc  test.  Data  are  average  ±  SD,  number  of  mice/experimental  group 
=  9-13- 

We  have  previously  shown  that  disabling  Nmpd  enhances  the  response  of  trabecular  bone  to  PTH  without 
compromising  the  cortical  compartments.  This  is  consistent  with  our  observation  that  Nmp4~/~  mice  from  the 
PTH 4- RAL  therpay  group  had  the  largest  cortical  bone  area  (Figures  4A  and  4B).  However,  there  was  no 
evidence  of  synergy  in  this  or  any  of  the  combination  treatments  for  this  parameter  nor  was  there  a  genotype  x 
treatment  effect. 

KEY  FINDING:  Nmp4~/~  mice  under  the  PTH+RAL  combination  therapy  harbored  more  bone  marrow 
osteoprogenitors  than  mice  under  the  other  treatments 


Percent  Bone  marrow  osteoprogenitors  (CD45-/CD105't‘/CD146**‘/nestin+) 
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Figure  5:  Percent  bone  marrow  (BM)  osteoprogenitors  as  evaluated  by 
CD457CDi05+/CDi46+/Nestin+  cells.  Statistical  analysis  used  a  non- 
parametric  Kruskal-Wallis  test  followed  by  a  Dunn’s  post  hoc.  See  text  for 
details. 


We  have  reported  that  healthy,  non- 
ovariectomized  Nmp4~/'  mice  exhibit  an 
increased  number  of  bone  marrow  (BM) 
osteoprogenitor  cells  irrespective  of  treatments 
We  also  demonstrated  that  the  percentage  of 
BM  CD457CDi05VCDi46VNestin+  cells 
paralleled  the  number  of  CFU-F®i'^p^°®+  cells  in 
our  miceS  During  the  first  year  of  this  DOD 
study  we  showed  that  this  enhanced  number  of 
Nmp4~/'  osteoprogenitors  wanes  I2wks  post¬ 
ovariectomy.  Here  we  report  that  this  increased 
pool  of  osteoprogenitors  is  maintained  i2wks 
post-op  in  null  mice  under  the  PTH -i- RAL 
therapy  (see  Figure  5).  A  Kruskal-Wallis 
analysis  with  Dunn’s  post  hoc  test  comparing 
the  %  BM  osteoprogenitors  in  the  WT-VEH 
cohort  with  all  the  other  groups  show  that  the 
Nmp4~/'  mice  under  the  PTH  4- RAL  therapy 
cohort  are  the  only  treatment  group  that  exhibit 
a  significant  difference  from  the  control. 
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KEY  FINDING:  Raloxifene  enhanced  WT  mesenchymal  stem/progenitor  cell  (MSPC)  mineralization 
We  have  shown  that  Nmp4~/~  expanded  MSPCs  exhibit  a  precocious  mineralization  (see  Appendix  Manuscript 
Number  ME-14-1406R1).  The  null  cells  begin  to  mineralize  by  Days  4-7  in  culture  whereas  the  WT  cells  have  a 
very  weak  mineralization  response.  There  are  a  number  of  studies  showing  that  osteoprogenitors  and 
osteoblasts  respond  directly  to  raloxifene^°'^5.  it  has  also  been  demonstrated  that  raloxifene  has  a 
mineralization-promoting  effect  on  female  mesenchymal  stem  cells^s.  As  part  of  Subtask  6.2:  Phenotype 
analyses  of  WT  and  Nmp4~/~  mesenchymal  stem  progenitor  cells  (MSPCs)  we  addressed  whether  raloxifene 
enhances  the  mineralization  response  to  ex  vivo  cultures  of  expanded  MSPCs.  Figure  6  outlines  the 
experimental  protocol  and  results 

These  cells  were  derived  from  bone  marrow  mononuclear  cells  of  male  mice  and  expanded  in  MesenCult  for  5 
passages  before  use  in  osteogenic  differentiation  studies.  Cells  are  used  between  passages  5  and  10^'’. 


y 


Day  0:  Seed  25K  cells  into  each  well  in  aMEM  complete  as  shown  in  the  Figure.  [03/16/15] 

Day  2:  Transfer  cells  to  osteogenic  medium  [VEH3]  or  RAL+osteogenic  medium. 

For  RAL  supplement,  add  1 .02pl  of  Raloxifene  Secondary  Stock  to  1  mis  of  complete  medium.  To  VEH3  add 
1.02pl  RAL  diluent/1  ml  of  complete  medium.  1584l*"  isisaR™  15841"^  isisrr^  15840^  isisrr"® 

Feed  the  cells  every  MWF.  '  ' 

Day  5,  Day  9,  and  Day  16,  stain  the  cells  with  alizarin  red. 


1584L'^T  1515RR'<° 

-RAL  +RAL  -RAL  +RAL 


1584LWT  1515RRKO 

-RAL  +RAL  -RAL  +RAL 


DAYS 


DAY  16 


Figure  6:  Raloxifene  (lOOnM)  clearly  enhances  WT  MSPC  mineralization.  Whether  it  has  an  impact  on  the  precocious  mineralization 
of  the  Nmp4~y  (KO)  remains  to  be  determined.  See  text  for  details 


We  made  the  following  observations  from  this  single  preliminary  experiment: 

•  At  Day  9,  post-seeding  the  Nmp4~/~  cells  show  a  strong  mineralization  response  in  both  the  presence 
and  absence  of  raloxifene  (lOonM),  but  the  WT  cells  have  little  to  no  evidence  of  mineralization. 

•  At  Day  16  the  WT  cells  in  raloxifene  are  showing  a  significant  increase  in  mineralization  compared  to 
the  cells  in  the  absence  of  this  drug 

•  The  Nmp4~/~  cells  have  mineralized  to  the  point  of  obscuring  any  potential  effect  of  the  drug  using  this 
qualitative  assay 
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ADDITIONAL  RESEARCH  ACTIVITIES; 

Breeding  of  conditional  knockout  (KO)  mice 

The  objectives  for  subtasks  4. 2. a  and  4.3. a  are  the  breeding  oiNmp^/A  3.6Col-Cre+,  Nmp4fl/fl  Cathepsin  K- 
Cre+,  and  Nmp4fl/fl-Cre-  mice.  After  breeding  our  germline  transmission  mice  with  the  B6.129S4- 
Gt(ROSA)26Sortum2(FLP)Sor/J  mice  to  remove  the  Neo  gene  we  obtained  evidence  that  the  5’  floxed  site  was 
absent  in  the  offspring.  We  are  currently  repeating  these  analyses  and  if  necessary  will  devise  a  strategy  for 
restoring  this  site. 

KEY  RESEARCH  ACCOMPLISHMENTS  FOR  2^0  YEAR: 

•  KEY  FINDING;  Mice  receiving  PTH+RAL  and  PTH+ZOL  showed  the  largest  BMD  increase  after  ovx. 

•  KEY  FINDING:  The  PTH+RAL  therapy  was  the  only  treatment  to  exhibit  a  synergistic  effect  on  total 
WB  and  spine  BMD 

•  KEY  FINDING:  Disabling  Nmp4  enhanced  PTH+RAL-induced  increases  in  femoral  BV/TV  and 
increased  the  synergistic  action  of  PTH  with  RAL  and  ZOL  in  both  the  femur  and  spine. 

•  KEY  FINDING:  Disabling  Nmp4  enhanced  the  distinction  between  the  PTH+RAL-induced  increase 
in  cortical  bone  area  and  the  other  combination  treatments 

•  KEY  FINDING:  Nmp4~/-  mice  under  the  PTH+RAL  combination  therapy  harbored  more  bone 
marrow  osteoprogenitors  than  mice  under  the  other  treatments 

•  KEY  FINDING:  Raloxifene  enhanced  WT  mesenchymal  stem/progenitor  cell  (MSPC)  mineralization 

REPORTABLE  OUTCOMES: 
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CONCLUSION :  Key  experimental  discovery  for  2°*^  year  of  DOD  study 

Our  key  discovery  is  that  disabling  Nmpy  significantly  improves  PTH+raloxifene  combination  therapy  in  an 
ovariectomized preclinical  osteoporosis  model  but  does  not  have  a  strong  impact  on  PTH+bisphosphonate 
combination  treatments 

“So  what?’ 

Women  comprise  the  fastest  growing  group  of  the  US  veterans  contributing  to  the  looming  osteoporosis 
epidemic  within  the  veteran  population.  The  Veterans  Affairs  (VA)  health  care  system  will  he  in  high  demand 
hy  female  veterans  of  Operation  Enduring  Freedom  and  Operation  Iraqi  Freedom.  PTH  is  the  only  FDA- 
approved  anabolic  osteoporosis  therapy  and  adds  significant  amounts  of  bone  to  the  osteoporotic  skeleton. 
Therefore,  this  drug  has  the  potential  to  restore  the  bone  lost  in  a  variety  of  VA  clinical  settings.  However,  a 
drawback  to  PTH  use  is  that  potency  declines  within  2  years  and  thus  it  is  not  suitable  as  a  long-term  therapy, 
which  is  problematic  in  treating  a  chronic  degenerative  disease.  Clinicians  have  attempted  to  improve  PTH 
therapy  by  adding  an  anti-catabolic  (e.g.  raloxifene  or  a  bisphosphonate)  to  the  treatment.  This  has  not  met 
with  success.  The  present  discovery  supports  our  contention  that  disabling  Nmpq  or  some  component  of  its 
pathway  will  unlock  the  block  on  PTH  combination  therapies  and  enhance/extend  regeneration  of  osteoporotic 
bone  in  post-menopausal  female  veterans. 
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Abstract: 

Parathyroid  hormone  (PTH)  is  an  osteoanabolic  for  treating  osteoporosis  but  its 
potency  wanes.  Disabling  the  transcription  factor  Nmp4  in  healthy,  ovary-intact  mice 
enhances  bone  response  to  PTH  and  BMP2  and  protects  from  unloading-induced 
osteopenia.  These  Nmp4-/-  mice  exhibit  expanded  bone  marrow  (BM)  populations  of 
osteoprogenitors  and  supporting  CD8+  T  cells.  To  determine  whether  the  Nmp4-/- 
phenotype  persists  in  an  osteoporosis  model  we  compared  PTH  response  in 
ovariectomized  (ovx)  wild  type  (WT)  and  Nmp4-/-  mice.  To  identify  potential  Nmp4 
target  genes  we  performed  bioinformatic/pathway  profiling  on  Nmp4  ChIP-seq  data. 
Mice  (12wks)  were  ovx  or  sham-operated  4wks  before  the  initiation  of  PTH  therapy. 
Skeletal  phenotype  analysis  included  pCT,  histomorphometry,  serum  profiles,  FACS 
sorting  and  the  growth/mineralization  of  cultured  WT  and  Nmp4-/-  BM  mesenchymal 
stem  progenitor  cells  (MSPCs).  ChIP-seq  data  were  derived  using  MC3T3-E1  pre¬ 
osteoblasts,  murine  embryonic  stem  cells,  and  two  blood  cell  lines.  Ovx  Nmp4-/-  mice 
exhibited  an  improved  response  to  PTH  coupled  with  elevated  numbers  of 
osteoprogenitors  and  CD8+  T  cells,  but  were  not  protected  from  ovx-induced  bone 
loss.  Cultured  Nmp4-/-  MSPCs  displayed  enhanced  proliferation  and  accelerated 
mineralization.  ChIP-seq/gene  ontology  analyses  identified  target  genes  likely  under 
Nmp4  control  as  enriched  for  negative  regulators  of  biosynthetic  processes. 

Interrogation  of  mRNA  transcripts  in  non-differentiating  and  osteogenic  differentiating 
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WT  and  Nmp4-/-  MSPCs  was  performed  on  90  Nmp4  target  genes  and  differentiation 
markers.  These  data  suggest  that  Nmp4  suppresses  bone  anabolism,  in  part,  by 
regulating  insulin-like  growth  factor  binding  protein  expression.  Changes  in  Nmp4 
status  may  iead  to  improvements  in  osteoprogenitor  response  to  therapeutic  cues. 
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Cover  Letter 


INDIANA  UNIVERSITY 


DEPARTMENT  OF  ANATOMY  AND  CELL  BIOLOGY 
Sclioul  ol  Medicine 


07/06/15 

Stephen  R.  Hammes,  M.D.,  Ph.D. 

Editor-in-Chief,  Molecular  Endocrinology 

Dear  Dr  Hammes; 

We  are  pleased  to  submit  our  revised  manuscript  Ref.:  Ms.  No.  ME-14-1406  “Genome-wide  mapping 
and  interrogation  of  the  Nmp4  anti-anabolic  bone  axis"  for  your  consideration.  We  thank  the  reviewers 
for  their  insightful  and  helpful  comments.  We  have  made  a  substantial  effort  to  comprehensively  address 
their  concerns  (see  separate  file  RESPONSE  TO  REVIEWERS).  This  includes  the  addition  of  a  study  to 
evaluate  gene  expression  in  non-differentiating  and  in  osteogenic-differentiating  WT  and  Nmp4'^^  MPSCs 
during  a  16-day  culture  period  at  five  time -points  using  a  custom  Taqman  Eow-Density  Array  system 
designed  for  96  senes  including  Nmp4  target  genes  identified  by  our  genome -wide  ChlP-seq  profiling, 
non-target  genes  that  drive  osteogenic  differentiation  and  marker  genes  of  this  process.  These  results 
provide  considerable  new  information  on  how  Nmp4  may  suppress  the  anabolic  response  to  PTH. 
Additionally,  these  new  data  strengthen  our  assertion  that  further  interrogation  of  the  Nmp4  anti-anabolic 
network  will  identify  pharmacologically  accessible  pathways  for  adding  new  bone  to  the  old  skeleton. 
The  supplemental  files  are  intended  for  online  publication.  Note  that  the  revised/new  text  is  in  red  font 
with  a  black  line  in  the  right  margin.  Finally,  there  is  no  requirement  for  a  PDB  number  and  we  are  not 
reporting  a  new  structure  or  compounds. 


Sincerely, 


Joseph  Bidwell 
Professor 


Medical  Science  Building  5035  635  Barnhill  Drive  Indianapolis,  IN  46202-5120  (317)  274-7494  tel  (317)  278-2040  fax 


Point-by-Point  Rebuttal 


RESPONSE  TO  REVIEWERS 


Ms.  No.  ME-14-1406  “Genome-wide  mapping  and  interrogation  of  the  Nmp4  anti-anabolic  bone  axis” 

We  thank  the  reviewers  for  the  helpful  and  insightful  comments.  We  believe  that  following  their  advice  has  made  this  a 
more  powerful  and  informative  study.  Please  note  that  the  revised/new  text  is  in  red  font  with  a  black  line  in  the  right 
margin. 

Comment:  The  phenot5^ic  studies  in  the  current  work  are  merely  confirmative 

Response:  The  reviewers  raise  an  important  issue  and  we  clearly  failed  to  emphasize  the  significance  of  the 
questions  addressed  by  these  experiments. 

•  Our  data  are  a  key  and  obligatory  preclinical  extension  of  the  Nmp4^^^  phenotype.  An  overriding  objective  of  the 
Nmp4  studies  is  to  determine  if  this  network  of  pathways  harbor  potential  targets  for  treating  postmenopausal 
osteoporosis.  Therefore  we  addressed  two  central  phenotype  questions  in  the  present  study: 

(i)  Does  disabling  Nmp4  protect  mice  from  ovariectomy -induced  bone  loss  as  it  does  in  disuse-induced 
osteopenia? 

(ii)  Does  the  exaggerated  response  to  PTH  therapy  in  Nmp4^^'  mice  persist  after  ovariectomy? 

Our  findings  that  the  null  mice  are  not  protected  from  ovarectomy-induced  bone  loss  yet  maintain  the  amplified 

response  to  PTH  therapy  is  an  important  discovery.  These  data  begin  to  demarcate  the  potential  Nmp4-based 
therapeutic  strategies.  We  have  revised  both  the  Introduction  and  Discussion  to  emphasize  these  points. 

Comment:  If  Nmp4  suppresses  Igfl  and  Bmp2,  and  cellular  synthetic  process,  why  there  is  no  difference 
between  WT  vs.  Nmp4-/-  mice  under  either  hasal  or  OVX-challenged  conditions?  Why  is  there  only  enhanced 
response  to  PTH  treatment? 

Response:  The  reviewer  raises  a  significant  question  about  features  of  the  Nmp4-/-  phenomenon.  Although  a 
definitive  answer  requires  further  study,  we  address  this  query  in  the  Discussion.  The  key  concepts  of  our 
response  follow: 

•  The  Nmp4  phenomenon  differs  from  the  results  of  recent  clinical  trials  in  which  neutralizing  sclerostin,  an 
inhibitor  of  the  Wnt  signaling  pathway  and  osteoblast  differentiation,  significantly  increases  baseline  bone 
mineral  density. 

•  Apparently  disabling  Nmp4  is  not  sufficient  for  driving  excess  bone  formation  but  instead  enlarges  the  number  of 
osteoprogenitors  that  can  be  recruited  once  activated  by  an  anabolic  cue.  The  downregulation  of  Cxclll  in  the 
Nmp4'^  cells  may  contribute  to  the  increased  number  of  CFU-F^”^p''°®^  osteoprogenitors  (see  low-density  array  data 
below). 

•  Once  activated  by  an  anabolic  cue,  e.g.  PTH,  the  Nmp4'''  cells  produce  autocrine/paracrine  factors  that  enhance 
the  replication  and  differentiation  of  neighboring  osteoprogenitors,  a  key  early  event  driving  the  PTH  anabolic 
response  (see  below). 

•  The  occasionally  observed  elevated  trabecular  volume  in  untreated  Nmp4'^^  mice  may  be  due  to  the  sporadic  local 
release  of  growth  factors,  e.g.,  Igfl  and/or  Bmp2. 

•  However  exogenous  pharmacological  doses  of  PTH  provide  a  stronger  stimulus  for  triggering  the  response 
leading  to  an  enhanced  bone  formation. 

•  The  enhanced  mineralization  in  Nmp4^'^  MSPC  culture  is  only  observed  upon  the  addition  of  the  anabolic  cue 
dexamethasone,  an  Nmp4  target  pathway. 


Comment:  Table  1,  why  didn’t  the  CTX  bone  resorption  marker  and  PINP  bone  formation  marker  show 
significant  differences  between  the  pre-op  and  post-op  group,  even  in  WT  mice? 


Response: 

•  Surgery-induced  changes  in  serum  bone  markers  took  longer  than  4wks  to  be  detected.  Our  mice  showed  the 
appropriate  ovx-induced  significant  increases  in  body  weight,  decreases  in  uterine  weight,  and  significant  bone 
loss  during  the  first  4wks  post-op.  Nevertheless,  our  study  was  not  a  longitudinal  design  in  which  repeated 
measures  on  each  animal  at  numerous  time  points  were  obtained,  thus  the  typical  small  increases  in  serum  CTX 
seen  in  other  ovx  mouse  studies  may  have  been  obscured.  A  more  extensive  time  course  for  harvesting 
histomorphometry  and  serum  samples  is  required  to  more  fully  characterize  the  anticipated  differences  in  WT  and 
Nmp4  '  dynamic  bone  remodeling.  However,  removal  of  a  single  outlier  in  the  serum  PINP  data  indicated  there 
was  a  21%  decline  in  serum  PINP  in  the  Nmp4''^  mice  during  the  first  4wks  post-op  but  no  significant  change  in 
the  WT  animals  (Genotype  p=0.37;  Treatment  p=0.16;  GxT  p=0.05,  2W  ANOVA,  see  Results  and  Table  1 
[revised]).  Additionally,  from  pre-op  to  post-op  12wks  (vehicle-treated  mice)  there  was  an  approximately  50% 
decline  in  PINP  (genotype  p=0.95;  treatment  p<0.0001;  GxT  p=0.36)  and  an  approximate  20%  decline  in  serum 
CTX  (genotype  p=0. 14;  treatment  p=0.0003;  GxT  p=0. 1 8).  This  suggests  a  decrease  in  bone  remodeling  in  the 
untreated  mice  predominantly  in  the  bone  formation  arm.  We  further  address  this  issue  in  the  revised  Discussion. 

Comment:  Figures  1  and  2,  The  Nmp4-/-  mice  without  PTH  seems  to  lose  more  bone  at  12wks  compared  to 
the  WT  without  PTH  treatment.  Why  does  this  happen? 

Comment:  The  difference  of  FEMUR  BV/TV  between  WT  and  Nmp4-/-  at  12  wks  seems  to  be  smaller  (or  the 
same)  comparing  to  the  difference  of  FEMUR  BV/TV  between  WT  and  Nmp4-/-  at  8 wks.  This  phenomenon 
also  can  be  found  for  the  L5  BV/TV.  Can  authors  explain  this? 

Response: 

•  We  agree  that  the  graphical  data  suggest  the  vehicle-treated  Nmp4^^^  mice  lose  more  bone  between  8-12wks  and 
that  the  difference  between  BV/TV  in  the  two  genotypes  is  smaller  at  the  later  time  point,  however  stringent 
statistical  analysis  does  not  support  this  interpretation.  Comparison  of  the  femoral  BV/TV  data  derived  from  the 
vehicle-treated  mice  at  the  8wk  and  12wk  time  points  indicates  a  significant  decline  in  both  genotypes  but  no 
genotype  or  genotype  x  time  interaction  (Genotype  p=0.44;  Time  p=  0.03;  GxT  0.06,  2W  ANOVA).  A  similar 
result  is  obtained  with  the  L5  BV/TV  data.  This  indicates  that  the  modest  enhanced  bone  loss  in  the  Amp4  ^  mice 
observed  4wks  post-op  was  stabilized  or  stabilizing  by  8wks  post-op. 

Comment:  I  wonder  if  the  difference  of  the  PTH  treatment  between  Nmp4-/-  and  WT  will  be  less  when  the 
follow-up  time  is  longer  than  12  weeks? 

Response: 

•  The  reviewer’s  query  whether  the  enhanced  PTH-induced  bone  formation  in  the  Nmp4^^'  mice  will  decline  after 
12wks  is  very  interesting.  Although  the  elevated  number  of  Amp4  ^  osteoprogenitors  declined  at  12wks  post-op, 
the  higher  number  of  Nmp4~''  CD8+  T  cells  did  not  decrease.  Therefore,  although  the  augmented  response  to  PTH 
might  be  weakened  in  these  mice,  it  may  not  disappear  since  the  persistent  increased  lymphocyte  number  would 
provide  extra  WntlOb.  We  have  revised  the  Discussion  to  incorporate  this  talking  point. 


Comment:  The  identification  of  Nmp4  target  genes  relied  solely  on  ChIP-seq  analysis.  As  binding  to  the 
promoter  may  lead  to  either  gene  activation  or  repression,  this  approach  can  only  reveal  the  binding 
properties,  but  does  not  reveal  functional  outcome.  Gene  expression  analysis  by  either  microarray  or  RNA- 
seq  in  parallel  would  be  more  informative. 

Comment:  The  findings  from  ChIP-seq  and  bioinformatics  analysis  should  be  more  extensively  validated  by 
ChIP-QPCR,  mutagenesis  and  functional  analysis,  followed  by  rescuing  experiments  to  determine  whether 
the  potential  targets  are  truly  functionally  significant. 

Comment:  Although  authors  spent  a  great  effort  to  create  chip-seq  data  and  the  data  is  very  interesting,  but, 
regarding  the  specific  genes/pathways  that  Nmp4  may  be  involved  in  PTH  treatment,  findings  from  chip- 
seq  data  seem  to  fall  short  and  is  inconclusive.  I  am  wondering  if  there  are  Nmp4  targeted  genes  that  are 
relevant  to  osteogenesis,  bone  mineralization,  PTH  and/or  PTH  treatment  relevant  pathways?  In  addition, 
have  authors  performed  transcriptome  profiling  in  relevant  cells/tissues  on  those  four  experiment  groups 
of  animals?  Did  IGFl  or  BMP  express  differently? 

Response: 

•  We  have  added  significantly  more  functional  data  to  the  revised  manuscript.  To  evaluate  gene  expression  in  non- 
differentiating  and  in  osteogenic-differentiating  WT  and  Nmp4^^^  MPSCs,  RNA  was  isolated  at  five  time  points 
from  both  genotypes  over  the  16-day  culture  period.  Individual  cDNAs  were  quantified  by  qRT-PCR  using  a 
custom  Taqman  Low-Density  Array  system  (Format  96a,  Applied  Biosystems,  Foster  City,  CA)  designed  for  96 
genes  including  Nmp4  target  genes  identified  by  our  genome -wide  ChlP-seq  profiling,  non-target  genes  that  drive 
osteogenic  differentiation  and  marker  genes  of  this  process.  Our  results  are  presented  in  Figures  9,  10  and 
Supplemental  Figure  2A-2E.  These  data  are  very  informative  and  suggest  a  mechanism  underlying  the  enhanced 
PTH-induced  anabolism  in  the  ^  mice.  Briefly,  the  Nmp4'^'  cells  showed  a  dramatic  surge  in  the  expression 
of  Igfbpl  and  downregulation  of  Igfbp4  during  the  early  stages  of  differentiation,  both  Nmp4  target  genes.  Igfbpl 
is  a  strong  autocrine/paracrine  anabolic  signal  and  a  recent  study  demonstrated  that  over-expression  of  this  gene 
accelerated  mineralization  and  enhanced  osteocalcin  (Bslap)  expression  in  bone  cells  [Xi  et  ah,  20141,  strikingly 
similar  to  the  phenotype  of  our  null  osteoprogenitors.  Igfbp4  is  a  potent  inhibitor  of  osteogenesis  [Miyakoshi  et  ah, 
19991.  Interestingly,  disabling  Nmp4  had  no  impact  on  the  mRNA  expression  of  the  target  genes  Igfl  and  Igfrl  but 
did  elevate  expression  of  Pdkl,  a  key  component  of  the  Igfl/insulin  signaling  pathway.  Near  the  end  of  the  16- 
day  differentiation  period  the  NmpP^~  cells  exhibited  an  enhanced  anabolic  profile  and  expressed  elevated  levels 
of  the  non-target  genes  Bmp2,  Pthlr,  and  Bslap.  We  also  described  the  markedly  downregulation  of  Cxclll  in  the 
Nmp4'^'  cells,  a  gene  that  plays  a  signiflcant  role  in  MSPC  osteolineage  commitment.  Other  differences  between 
the  WT  and  NmpP'^  mRNA  expression  profiles  are  described.  We  agree  that  mutagenesis  and  rescuing 
experiments  are  the  next  step  in  analyzing  the  Nmp4  network  but  we  respectfully  contend  that  they  are  beyond 
the  scope  of  the  present  study. 

•  Our  Nmp4  ChlP-seq-derived  map  is  the  vital  first  step  for  deciphering  the  gene  regulatory  networks  underlying 
the  biological  processes  of  this  novel  anti-anabolic  axis.  Our  low-density  array  results  are  the  first  step  in 
validating  the  reliability  of  this  map. 


Comment:  Tables  1  and  2:  It  would  be  great  to  know  the  number  of  animals  in  each  group. 
Response:  We  have  placed  this  information  in  Tables  1  and  2. 
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ABSTRACT: 


Parathyroid  hormone  (PTH)  is  an  osteoanabolic  for  treating  osteoporosis  but  its  potency  wanes.  Disabling 
the  transcription  factor  Nmp4  in  healthy,  ovary-intact  mice  enhances  bone  response  to  PTH  and  BMP2 
and  protects  from  unloading-induced  osteopenia.  These  Nmp4~^'  mice  exhibit  expanded  bone  marrow 
(BM)  populations  of  osteoprogenitors  and  supporting  CD8^  T  cells.  To  determine  whether  the  Nmp4'^^ 
phenotype  persists  in  an  osteoporosis  model  we  compared  PTH  response  in  ovariectomized  (ovx)  wild 
type  (WT)  and  Nmp4''^  mice.  To  identify  potential  Nmp4  target  genes  we  performed 
bioinformatic/pathway  profiling  on  Nmp4  ChIP-seq  data.  Mice  (12wks)  were  ovx  or  sham-operated  4wks 
before  the  initiation  of  PTH  therapy.  Skeletal  phenotype  analysis  included  pCT,  histomorphometry, 
serum  profiles,  FACS  sorting  and  the  growth/mineralization  of  cultured  WT  and  Am/74  ^  BM 
mesenchymal  stem  progenitor  cells  (MSPCs).  ChIP-seq  data  were  derived  using  MC3T3-E1  pre¬ 
osteoblasts,  murine  embryonic  stem  cells,  and  two  blood  cell  lines.  Ovx  Am/74  mice  exhibited  an 
improved  response  to  PTH  coupled  with  elevated  numbers  of  osteoprogenitors  and  CD8^  T  cells,  but 
were  not  protected  from  ovx -induced  bone  loss.  Cultured  Am/74  ^  MSPCs  displayed  enhanced 
proliferation  and  accelerated  mineralization.  ChIP-seq/gene  ontology  analyses  identified  target  genes 
likely  under  Nmp4  control  as  enriched  for  negative  regulators  of  biosynthetic  processes.  Interrogation  of 
ruRNA  transcripts  in  non-differentiating  and  osteogenic  differentiating  WT  and  Nmp4'^^  MSPCs  was 
performed  on  90  Nmp4  target  genes  and  differentiation  markers.  These  data  suggest  that  Nmp4 
suppresses  bone  anabolism,  in  part,  by  regulating  insulin-like  growth  factor  binding  protein  expression. 
Changes  in  Nmp4  status  may  lead  to  improvements  in  osteoprogenitor  response  to  therapeutic  cues. 

250  words 
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64 


INTRODUCTION: 


65  Patients  with  severe  osteoporosis  are  often  treated  with  parathyroid  hormone  (PTH),  a  potent 

66  osteoanabolic  agent  [Kraenzlin  &  Meier,  2011],  however,  the  bone-building  ability  of  this  drug  or  its 

67  ‘anabolic  window’  wanes,  likely  due  to  latent  increases  in  bone  resorption.  [Cipriani  et  ah,  2012;  Baron  & 

68  Hesse,  2012;  Yu  et  ah,  2011],  This  limits  its  effectiveness  to  treat  a  chronic  degenerative  disease.  Recent 

69  advances  in  bone-forming  agents  have  shown  that  one  can  increase  the  extent  of  bone  mass  accrual  with 

70  anti-SOST  treatment  compared  to  PTH  (McClung  et  ah,  2014).  However,  there  may  be  unique  pathways 

7 1  triggered  by  PTH,  which  allows  for  sustained  targeting  of  early  osteogenesis  as  evidenced  by  serum 

72  markers  of  bone  formation  such  as  N-terminal  propeptide  of  type  1  procollagen  (PINP)  and  osteocalcin 

73  (OCN,  Padhi  et  ah,  2011;  Saag  et  ah,  2009).  In  contrast  to  PTH,  anti-SOST  antibodies  may  have  a  limited 

74  capacity  for  targeting  osteoprogenitors  as  evidenced  by  a  relatively  transient  up-regulation  of  collagen- 

75  based  markers  such  as  PINP  (McClung  et  ah,  2014).  Therefore  given  PTH’s  unique  mode  of  action, 

7 6  therapies  that  could  enhance  PTH-mediated  recruitment  of  osteoprogenitors  may  add  value  to  some 

77  patients.  How  to  achieve  this  enhancement  is  not  clear.  For  example,  attempts  to  extend  and  enhance  PTH 

78  efficacy  by  combining  treatment  with  anti-resorptive  medications  have  met  with  mixed  success  and  have 

79  generally  been  underwhelming  [Cosman  et  ah,  2011;  Finkelstein  et  ah,  2010;  Black  et  ah,  2003]. 

80  Blocking  the  activity  of  Nmp4/ClZ  (nuclear  matrix  protein  4/cas  interacting  zinc  finger  protein, 

81  ‘Nmp4’)  in  mice  dramatically  enhanced  their  response  to  anabolic  doses  of  PTH  [He  et  ah,  2013; 

82  Childress  et  ah,  2011;  Robling  et  ah,  2009],  suggesting  a  potential  strategy  for  an  adjuvant  therapy  [Krane, 

83  2005].  Intermittent  exogenous  doses  of  hormone  stimulated  equivalent  new  bone  formation  in  wild  type 

84  (WT)  and  ^  mice  during  the  first  2wks  of  challenge,  but  at  3wks  of  treatment  the  null  mice 

85  exhibited  greater  than  a  2-fold  increase  in  new  trabecular  bone  compared  to  their  WT  littermates 

86  [Childress  et  ah,  2011].  This  augmented  skeletogenesis  in  the  Nmp4^^^  mice  was  extended  to  7wks  of 

87  treatment  and  was  observed  in  the  femur,  tibia,  and  vertebra.  Serum  osteocalcin  continued  to  rise  at  this 

88  time  point  in  the  ^  mice  but  had  decreased  in  the  WT  animals  [Childress  et  ah,  2011].  However,  the 

89  PTH  response  of  the  cortical  compartment  was  equivalent  throughout  treatment  in  the  WT  and  null  mice 
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[Robling  et  al.,  2009].  This  suggests  that  disabling  Nmp4  accelerates  and  enhances  the  response  of  bone 
to  intermittent  PTH  [Childress  et  ah,  2011]. 

Nmp4''^  bone  may  have  a  generalized  accelerated  and  heightened  response  to  systemic  or  local 
anabolic  cues.  For  example,  these  mice  also  exhibited  augmented  BMP2-induced  ectopic  bone  formation 
compared  to  their  WT  littermates  [Morinobu  et  al.,  2005].  The  Nmp4-Wi\\  mice  showed  an  accelerated 
osseous  regeneration  after  marrow  ablation  [Morinobu  et  al.,  2005]  and  did  not  lose  bone  during  hind 
limb  unloading,  which  appeared  to  derive  from  an  enhanced  osteoblast  activity  [Hino  et  al.,  2007]. 

Prerequisite  for  an  adjuvant  therapy  target,  disabling  Nmp4  has  little  impact  on  the  health, 
longevity,  or  global  baseline  phenotype  of  the  mouse,  with  a  few  exceptions.  The  Nmp4^^^  baseline 
skeletal  phenotype  (i.e.,  bone  mineral  density  and/or  content  and  trabecular  architecture)  is  generally 
equivalent  compared  to  WT  animals],  although  we  have  occasionally  observed  an  unprovoked  increase  in 
bone  properties  in  Nmp4'^^  mice  [He  et  al.,  2013,  Childress  et  al.,  2011;  Robling  et  al.,  2009;  Morinobu  et 
al.,  2005].  Similarly,  male  Nmp4^^^  mice  exhibit  variable  degrees  of  spermatogenic  cell  degeneration 
resembling  germinal-cell  aplasia  with  focal  spermatogenesis  resulting  in  sporadic  infertility  [Nakamoto  et 
al,  2004]. 

Our  recent  work  suggests  that  the  cellular  basis  of  the  osteoanabolic  repressor  function  of  Nmp4 
is  due  to  its  effect  on  the  bone  marrow  derived  stromal  stem/progenitor  cells  aka  mesenchymal  stem 
progenitor  cells  (MSPCs)  [He  et  al.,  2013].  Nmp4''^  mice  have  significantly  more  osteoprogenitor  cells  in 
their  marrow,  which  lie  in  wait  to  be  quickly  mobilized  to  differentiate  into  active  osteoblasts  upon 
stimulation  with  various  osteoanabolic  stimuli  [He  et  al.,  2013].  There  was  no  difference  between  WT  and 
Nmp4  '^  BM  cellularity  or  profiles  of  several  blood  elements  however,  the  null  mouse  exhibited  a  4-fold 
increase  in  CD457CD105VnestinVCD146^  BM  osteoprogenitor  cells.  These  markers  are  a  common 
hallmark  to  CFU-F  cells  with  osteogenic  potential  [Isem  et  al.,  2013;  Mendez-Ferrer  et  al.,  2010]  and 
indeed  4-fold  more  CFU-F^''^p’'“^  and  CFU-F°'’  cells  have  been  recovered  from  these  mice  compared  to 
the  WT  animals  [He  et  al.,  2013;  Morinobu  et  al.,  2005].  A  second,  related  phenomenon  we  have 
observed  in  Nmp4^^'  mice  is  a  2-fold  increase  in  the  prevalence  of  CD8^  T-cells  in  the  femoral  marrow — 
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the  lymphocyte  population  that  provides  potent  input  to  induce  MSPCs  down  the  osteoblast 
differentiation  pathway  [He  et  ah,  2013;  Bedi  et  ah,  2012;  Terauchi  et  ah,  2009;  Li  et  ah,  2014].  These 
blood  cells  express  the  PTHRl  receptor  and  support  the  PTH  anabolic  response  via  the  release  of  WntlOb 
upon  hormone  challenge,  which  drives  osteoprogenitor  differentiation  to  pre -osteoblasts  and  mature 
matrix -producing  bone  cells  [Bedi  et  ah,  2012;  Terauchi  et  ah,  2009;  Li  et  al.,  2014], 

There  is  little  information  on  the  molecular  mechanisms  and  cellular  pathways  that  mediate  the 
anti-anabolic  action  of  Nmp4.  This  transcription  factor  is  a  Cys2His2  zinc  finger  protein  that  primarily 
localizes  to  the  nucleus  although  there  is  evidence  for  cytoplasmic  activity  [Bidwell  et  al.,  2012; 
Nakamoto  et  al.,  2000].  The  zinc  fingers  recognize  the  DNA  minor  groove  of  an  AT -rich  consensus 
sequence  and  two  transactivation  domains  can  suppress  or  activate  transcription  depending  on  the  cellular 
context  [Shah  et  al.,  2004;  Torrungruang  et  al.,  2002;  Thunyakitpisal  et  al.,  2001;  Nakamoto  et  al.,  2000; 
Alvarez  et  al.,  1998].  The  amino  terminus  of  the  rodent  protein  contains  an  SH3-binding  domain  that 
associates  with  the  adaptor  signaling  protein  pi 30Cas  [Nakaomoto  et  al.,  2000],  but  the  functional 
significance  of  this  interaction  remains  unknown. 

The  Nmp4'''  progenitor  cells  and  their  progeny  have  an  exaggerated  stimulus  response  at  the 
levels  of  transcription  and  cell  signaling  [Alvarez  et  al.,  2012;  Yang  et  al.,  2010;  Shen  et  al.,  2002]. 
Nmp4-Wi\\  bone  marrow  stromal  cells  (BMSCs)  show  an  enhanced  transcriptional  response  to  PTH  and 
BMP2  [Alvarez  et  al.,  2012;  Shah  et  al.,  2004;  Shen  et  al.,  2002].  The  Nmp4^'^  derived  calvarial  cells 
exhibit  an  increased  load-induced  phosphorylation  of  Pi3k  and  Akt  and  beta-catenin  nuclear  translocation 
[Yang  et  al.,  2010].  Analogous  to  heightened  response  to  anabolic  signals  mNmp4^^'  osteolineage  cells, 
osteoclast  preparations  from  the  null  mice  exhibited  a  heightened  response  to  the  remodeling  signals  of 
RANKL  and  M-CSF  (Childress  et  al.,  2011). 

Two  essential  genotype-phenotype  questions  remaining  to  be  addressed  are  (i)  whether  the  Nmp4- 
null  mouse  is  resistant  to  ovariectomy  (ovx)-induced  bone  loss  and  (ii)  if  disabling  Nmp4  improves  PTH- 
based  bone  therapy  in  the  OVX  model.  This  is  a  focal  preclinical  extension  of  the  Am/74  ^  phenotype 
necessary  before  this  gene  and  its  associated  pathways  can  be  considered  potential  targets  for  an  adjuvant 
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therapy.  Additionally,  we  used  expanded  cultures  of  WT  and  Nmp4'^  mesenchymal  stem/progenitor  cells 
(MSPCs)  to  probe  the  cell  autonomous  proliferative  and  mineralization  activities  of  this  cell  population. 
To  delineate  the  framework  of  the  Nmp4  anti-anabolic  network  we  performed  genome -wide  chromatin 
immunoprecipitation  sequencing  (ChIP-seq)  on  MC3T3-E1  cells  and  combined  these  data  with  the  data 
available  for  Nmp4  (a.k.a.  Znf384)  from  the  Mouse  Encyclopedia  of  DNA  Elements  (ENCODE) 
Consortium  for  transcription  factors  [Stamatoyannopoulos  et  ah,  2012].  Bioinformatic  profiling,  gene 
ontology  (GO),  and  pathway  analysis  were  performed  on  these  data  sets  to  infer  a  map  of  the  negative 
regulation  of  bone  anabolism  under  Nmp4  control.  Interrogation  of  mRNA  transcripts  in  non- 
differentiating  and  osteogenic  differentiating  WT  and  Nmp4^^^  MSPCs  was  performed  on  90  Nmp4  target 
genes  and  differentiation  markers  to  inaugurate  validation  of  the  Nmp4  anti-anabolic  network. 

MATERIALS  AND  METHODS: 

Mice:  Male  and  female  Nmp4^^^  mice,  backcrossed  onto  a  C57BL/6J  background  for  7  generations  [He  et 
ah,  2013,  Childress  et  ah,  2011;  Robling  et  ah,  2009],  and  their  WT  littermates  were  produced  and 
maintained  in  our  colony  at  Indiana  University  Bioresearch  Facility,  Indiana  University  School  of 
Dentistry.  Our  local  Institutional  Animal  Care  and  Use  Committee  approved  all  husbandry  practices  and 
experimental  procedures  and  regimens  described  in  this  investigation. 

Bilateral  ovariectomy  surgery:  12wk-old  virgin  mice  were  anesthetized  using  isoflurane  inhalation 
followed  by  a  mixture  of  xylazine  and  ketamine  administered  intraperitoneally.  A  1  -2cm  dorsal  incision 
was  made  in  the  midline  below  the  level  of  the  last  rib  and  the  skin  bluntly  dissected  from  the  muscle  on 
either  side  of  the  incision.  Through  the  skin  incision,  the  muscle  wall  was  incised  1  cm  lateral  to  the 
midline  1  -2cm  below  the  last  rib  to  enter  the  abdominal  cavity.  The  periovarian  fat  pad  was  located  and 
gently  grasped  and  exteriorized.  Care  was  taken  not  to  directly  handle  the  ovary  to  avoid  abdominal 
implantation  of  ovarian  tissue.  While  holding  the  periovarian  fat  pad  with  forceps,  the  fallopian  tube 
between  the  fat  pad  and  uterus  was  clamped  and  crushed  using  mosquito  hemostats.  The  crushed  area  was 
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168  cut  with  scissors  and  the  fat  pad  with  ovary  removed.  The  procedure  was  repeated  on  the  contralateral 

169  side.  The  skin  incision  was  closed  with  one  or  two  surgical  wound  clips.  The  sham  surgeries  involved  all 

170  the  outlined  steps  except  the  crushing  the  fallopian  tubes  and  the  actual  removal  of  the  ovaries.  To 

171  confirm  the  efficacy  of  OVX,  uteri  were  weighed  following  euthanasia. 

172 

173  PTH  treatment:  At  16  wks  of  age,  ovx  animals  were  sorted  into  four  treatment  groups  based  on  equivalent 

174  mean-group-body  weight.  These  four  groups  included  1)  vehicle-treated  WT;  2)  PTH-treated  WT;  3) 

175  vehicle-treated  NmpP^^  and  4)  PTH-treated  Nmp4  '^  mice.  Mice  were  injected  subcutaneously  (sc)  with 

176  synthetic  human  PTH  1-34  acetate  salt  (Bachem  Bioscience  Inc,  PA)  at  30pg/kg/day,  daily  or  vehicle 

177  control  (0.2%  BSA/l.OpN  HCl  in  saline,  Abbott  Laboratory,  North  Chicago,  IL)  for  the  length  of  time 

178  indicated. 

179 

180  Cell  culture:  Cells  from  ATCC  (MC3T3-E1  subclone  4)  were  maintained  in  a-MEM  medium 

181  supplemented  with  100  lU/ml  penicillin,  100  pg/ml  streptomycin,  25  pg/ml  amphotericin,  2  mM  L- 

182  glutamine  (Gibco  BRL,  Grand  Island,  NY),  ascorbic  acid  (50pg/ml,  Sigma-Aldrich,  St  Louis,  MO),  and 

183  10%  fetal  bovine  serum  (FBS;  Sigma-Aldrich,  St  Louis,  MO).  Expanded  mesenchymal  stem/progenitor 

184  cell  (MSPC)  cultures  were  established  as  previously  described  [Wu  et  ah,  2006].  Briefly,  long  bone  BM 

185  was  isolated  from  euthanized  mice  6-8wks  of  age  and  the  mononuclear  cells  (BMMNCs)  were  isolated 

186  using  a  Ficoll  gradient.  These  cells  were  plated  in  Mesencult™  Media  -i-  Mesencult™  Stimulatory 

187  Supplement  (StemCell™  Technologies,  Vancouver  BC,  Canada)  and  maintained  in  culture  for  3-4wks 

188  without  passage  and  fed  every  5-7  days  by  removing  50%  of  the  old  media  and  adding  50%  fresh  media, 

189  very  gently  so  as  not  to  disturb  the  cells.  At  approximately  80%  confluence,  the  cells  were  passaged  at  1 :3 

190  dilution  for  two  more  passages  before  use  or  were  frozen  for  storage.  Cells  were  used  for  experiments 

191  between  passages  5-10.  For  comparing  cell  proliferation  rates  between  WT  and  NmpP^  MSPCs,  the  cells 

192  were  transferred  to  a-MEM  medium  without  the  ascorbic  acid  in  12-well  plates  at  5,000  cells/well  (Day 
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0).  Cells  were  counted  on  Day  2,  4,  and  6  post-seeding  prior  to  refreshing  the  medium  for  the  remaining 
cells.  To  evaluate  mineralizing  capacity  cells  were  transferred  to  osteogenic  differentiation  medium  and 
after  48hrs  (Day  0),  which  was  comprised  of  a-MEM  supplemented  with  ascorbic  acid  (5-50pg/ml, 

Sigma  Aldrich),  dexamethasone  (0-1  OnM,  Sigma-Aldrich),  and  lOmM  glycerol  2-phosphate  disodium 
salt  hydrate  (BGP,  Sigma-Aldrich).  For  controls,  cells  were  passaged  into  fresh  MesenCult  medium 
without  the  osteogenic/mineralization  supplements.  Cells  were  stained  for  alkaline  phosphatase  activity 
using  naphthol  AS-MX  phosphate  and  fast  red  violet  B  salt  following  the  manufacturer's  instructions 
(Sigma  cat#  85L3R-1KT)  or  for  mineralization  using  alizarin  red. 

To  compare  mRNA  expression  profiles  of  select  genes  in  non-differentiating  and  osteogenic- 
differentiating  WT  and  Nmp4''  MSPCs,  cells  were  seeded  into  12-well  plates  at  either  10,000  or  25,000 
cells/well  in  Mesencult™  Media  +  Mesencult™  Stimulatory  Supplement.  Those  cells  at  the  lower  seeding 
density  were  harvested  on  Day  3  post-seeding  (non-differentiating).  The  remaining  cells  were  transferred 
to  osteogenic  differentiation  medium  48hrs  post-seeding  and  harvested  on  Days  5,  7,  9,  and  16. 

Flow  cytometry.  Cellular  surface  marker  profiles  from  BM  and  peripheral  blood  (PBL)  were  assessed  as 
previously  described  [He  et  ah,  2013].  The  antibodies  employed  for  flow  cytometry  were  obtained  from 
BD  Biosciences  (San  Jose,  CA).  Stained  cells  were  analyzed  on  an  FACS  Calibur  (BD  Biosciences)  and 
results  were  quantified  using  FlowJo  Version  8.8.6  software  (TreeStar  Inc,  Ashland  OR). 

Micro  computed  tomography  (pCT):  Trabecular  bone  architecture  was  analyzed  as  we  have  previously 
described  [He  et  ah,  2013;  Childress  et  ah,  2011].  Briefly,  femurs  and  F5  vertebra  were  excised  from  the 
WT  and  Nmp4  '^  mice  after  euthanasia,  the  muscle  and  connective  tissue  removed,  and  the  bones 
transferred  to  10%  buffered  formalin,  4°C  for  48  hr,  after  which  the  bones  were  placed  in  70%  ethanol 
(4°C)  until  analyzed.  For  femur  analysis  a  2.6-mm  span  (~5  mm^  of  medullary  space)  of  the  excised  distal 
femoral  metaphysis  was  scanned  in  70%  ethanol  on  a  desktop  pCT  (pCT  35;  Scanco  Medical  AG, 
Bassersdorf,  Switzerland)  at  10  pm  resolution  using  55-kVp  tube  potential  and  400-msec  integration  time. 
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219  to  measure  three-dimensional  morphometric  properties.  The  entire  vertebra  (L5)  were  scanned  using 

220  standard  methods  (Skyscan  1172).  Bones  were  reconstructed  and  analyzed  using  the  manufacturer’s 

221  software.  The  trabecular  bone  between  the  two  growth  plates  was  isolated  from  the  cortical  shell  via 

222  manual  tracing  and  assessed  for  trabecular  architecture.  From  the  three  dimensional  reconstructions  the 

223  following  parameters  were  obtained  using  the  Scanco  and  Skyscan  software  analyses:  trabecular  bone 

224  volume  per  total  volume  (BV/TV,  %),  connectivity  density  (Conn.D,  mm'^),  structure  model  index  (SMI), 

225  trabecular  number  (Tb.N,  mm"'),  trabecular  thickness  (Tb.Th,  mm),  and  spacing  (Tb.Sp,  mm)  [Bouxsein 

226  etal.,  2010]. 

227 

228  Bone  histomorphometry:  All  histomorphometric  parameters  were  obtained  as  previously  described 

229  [Childress  et  ah,  2011]  following  the  ASBMR  guidelines  [Dempster  et  ah,  2013].  Briefly,  mice  were 

230  administered  intraperitoneal  injections  of  calcein  green  (20  mg/kg;  Sigma -Aldrich)  and  alizarin  red  (25 

231  mg/kg,  Sigma-Aldrich)  6  and  3  days  before  euthanasia,  respectively.  The  femur  marrow  cavity  was 

232  exposed  via  cutting  the  anterior  face  of  the  epiphyseal  plate.  Bones  were  embedded  in  methyl- 

233  methacrylate  subsequent  to  dehydration  with  graded  alcohols,  sectioned  (4pm)  with  a  Leica  RM2255 

234  microtome  (Leica  Microsystems,  Wetzlar,  Germany),  and  mounted  unstained  on  microscope  slides  and 

235  imaged  under  fluorescent  light  with  a  microscope  system  [Childress  et  ah,  2011].  Bone  formation  rate 

236  (BFR),  mineral  apposition  rate  (MAR),  and  mineralizing  surface  (MS/BS)  were  obtained  from  a  0.03mm^ 

237  metaphyseal  region  of  interest  from  250pm  to  1750pm  below  the  growth  plate  using  ImagePro  3.1 

238  software  (Media  Cybernetics,  Bethesda,  MD,  USA). 

239 

240  Serum  biochemistry.  We  analyzed  serum  N-terminal  propeptide  of  type  1  procollagen  (PINP)  to  evaluate 

241  global  bone  formation  in  our  experimental  mice  using  the  Rat/Mouse  PINP  ElA  from  IDS 

242  Immunodiagnostic  Systems  (Scottsdale,  AZ)  following  the  manufacturer’s  instructions.  To  follow  bone 

243  resorption  we  analyzed  serum  C-terminal  telopeptides  (CTX)  with  the  RatLaps^"^  ELISA 

244  (Immunodiagnostic  Systems  Inc)  [Childress  et  ah,  2011]. 
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245 

246  Quantitative  real-time  PCR  (qRT-PCR)  analysis:  ChIP-qPCR  was  used  to  authenticate  select  ChlP-seq 

247  profiles  employing  SYBR  Green  assays  and  SYBR  Green  Supermix  (Bio-rad,  Hercules,  CA).  qRT-PCR 

248  reactions  were  carried  out  in  triplicate  on  specific  genomic  regions.  The  resulting  signals  were  normalized 

249  for  primer  efficiency  by  carrying  out  qRT-PCR  reactions  for  each  primer  pair  using  Input  DNA. 

250  To  evaluate  gene  expression  in  non-differentiating  and  in  osteogenic-differentiating  WT  and  Nmp4'^' 

251  MPSCs,  RNA  was  isolated  with  RNAeasy  columns  according  to  the  manufacturer’s  instructions  (Qiagen, 

252  Gaithersburg,  MD).  The  RNA  was  reverse-transcribed  via  the  High-Capacity  cDNA  Reverse 

253  Transcription  Kit  (Applied  Biosystems,  Foster  City,  CA).  RNA  expression  profiling  was  performed  on 

254  three-four  replicates  per  time  point  for  both  genotypes  over  the  16-day  culture  period.  Individual  cDNAs 

255  were  quantified  by  qRT-PCR  using  a  custom  TLDA  system  (Format  96a,  Applied  Biosystems,  Foster 

256  City,  CA)  designed  for  96  genes  including  Nmp4  target  genes  identified  by  our  genome -wide  ChIP-seq 

257  profiling,  osteogenic  differentiation  marker  genes,  and  candidate  normalizer  genes.  All  experiments  were 

258  performed  in  biological  quadruplicate  or  triplicate  with  TaqMan  fast  advanced  master  mix  (Applied 

259  Biosystems)  on  a  QuantStudio™  7  Flex  Real-Time  PCR  System.  The  probes  used  are  listed  in 

260  Supplemental  Table  1.  We  used  the  ExpressionSuite  vl.0.4™  analysis  software  (Applied  Biosystems)  to 

261  analyze  these  data.  This  software  utilizes  the  comparative  threshold  cycle  (A ACT)  method  to  quantify 

262  relative  gene  expression  across  a  large  number  of  genes  and  samples.  The  software  provides  options  to 

263  normalize  expression  data  using  either  global  normalization  or  endogenous  controls  and  calculates  fold 

264  changes  with  P  values.  Gene  expression  data  were  normalized  to  five  endogenous  controls  (18S,  Gusb, 

265  Rplp2,  B2m,  and  Hprt)  although  we  report  Gusb  and  Rplp2  data  here.  In  all  experiments,  the  CT  upper 

266  limit  was  set  to  40,  meaning  that  all  mRNA  detectors  with  a  CT  value  greater  than  or  equal  to  40  were 

267  excluded.  The  multiple-comparisons  correction  (Benjamini-Hochberg  method  for  false-discovery  rate) 

268  was  applied  to  the  data  and  a  P  value  of  <0.05  was  considered  significant.  Additionally,  individual  qRT- 

269  PCR  reactions  were  performed  to  monitor  the  expressions  of  Sp7  (osterix,  Mm00504574_ml)  and  Bglap 

270  (osteocalcin,  Mm03413826_mH)  using  Rplp2  as  the  normalizer  (Mm03059047_gH).  The  prepared 
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cDNA  was  used  to  set  up  qRT-PCR  reactions  using  FastStart  Universal  Probe  Master  mix  (Rox)  (Roche 


Life  Science,  Indianapolis,  IN). 

Chromatin  immunoprecipitation  sequencing  (ChIP-seq)  and  CHIP  analysis:  Cells  from  ATCC  (MC3T3- 
E1  subclone  4)  were  seeded  into  twenty-one  150mm  plates  at  an  initial  density  of  50,000  cells/plate  (320 
cells/cm^)  and  maintained  in  aMEM  complete  medium  -l-  ascorbic  acid.  On  Day  14  post-seeding,  cells 
were  treated  with  25nM  hPTH(l-34)  or  vehicle  control  for  Ihr  before  harvest.  Subsequent  to  treatment 
cells  were  fixed  with  1%  formaldehyde  for  1 5min  and  quenched  with  0.125M  glycine.  Cell  pellets  were 
frozen  in  an  ethanol  dry  ice  bath  and  shipped  to  Active  Motif  for  FactorPath'^'^  analysis.  The  chromatin 
was  isolated  from  the  pellets  by  adding  lysis  buffer  followed  by  disruption  with  a  Bounce  homogenizer. 
Eysates  were  sonicated  and  the  DNA  sheared  to  an  average  length  of  300-500  bp.  Genomic  DNA  (Input) 
was  prepared  by  treating  aliquots  of  chromatin  with  RNase,  proteinase  K  and  heat  for  de -crosslinking, 
followed  by  ethanol  precipitation.  Pellets  were  resuspended  and  the  resulting  DNA  was  quantified  on  a 
NanoDrop  spectrophotometer.  Extrapolation  to  the  original  chromatin  volume  allowed  quantitation  of  the 
total  chromatin  yield.  An  aliquot  of  chromatin  (30pg)  was  precleared  with  protein  A  agarose  beads 
(Invitrogen,  ThermoFisher  Scientific,  Waltham,  MA).  Genomic  DNA  regions  of  interest  were  isolated 
using  4pg  antibody  against  ZNF384  (Sigma  FtPA004051,  Eot  A57874).  Complexes  were  washed,  eluted 
from  the  beads  with  SDS  buffer,  and  subjected  to  RNase  and  proteinase  K  treatment.  Crosslinks  were 
reversed  by  incubation  overnight  at  65  °C,  and  ChIP  DNA  was  purified  by  phenol-chloroform  extraction 
and  ethanol  precipitation. 

Chip  Sequencing  (Illumina):  ChIP  and  Input  DNAs  were  prepared  for  amplification  by  converting 
overhangs  into  phosphorylated  blunt  ends  and  adding  an  adenine  to  the  3 ’-ends.  Illumina  genomic 
adapters  were  ligated  and  the  sample  was  size-fractionated  (200-300  bp)  on  an  agarose  gel.  After  a  final 
PCR  amplification  step  (18  cycles),  the  resulting  DNA  libraries  were  quantified  and  sequenced  on  HiSeq 
2000.  Sequences  (50nt  reads,  single  end)  were  aligned  to  the  mouse  genome  (mmlO)  using  the  BWA 


12 


297  algorithm.  Alignments  were  extended  in  silieo  at  their  3 ’-ends  to  a  length  of  150  bp,  which  is  the  average 

298  genomic  fragment  length  in  the  size-selected  library,  and  assigned  to  32-nt  bins  along  the  genome.  The 

299  resulting  histograms  (genomic  “signal  maps”)  were  stored  in  BAR  and  bigWig  files.  ZNF384  peak 

300  locations  were  determined  using  the  MACS  algorithm  (vl.4.2)  with  a  cutoff  of  pvalue  =  le-7  [Li  and 

301  Durbin,  2009]. 

302 

303  Bioinformatic  profiling:  In  addition  to  generating  our  own  Nmp4  ChIP-seq  data  from  the  MC3T3-E1  cells 

304  we  used  Nmp4  (Znf384)  ChlP-seq  data  from  murine  embryonic  stem  cell  line  (ES-E14)  and  the  B-cell 

305  lymphoma  cell  lines  Chl2  and  MEE  from  the  ENCODE  Consortium  for  transcription  factors  2011  Freeze 

306  data  sets  in  NarrowPeak  format  (Rosenbloom  et  ah,  2013).  To  assign  an  Nmp4  peak  to  a  promoter  region 

307  it  had  to  be  within  -5kb  to  -l-2kb  from  a  transcription  start  site  (TSS).  To  assign  a  peak  to  an  intragenic 

308  region  it  had  to  be  located  within  the  range  defined  by  the  TSS  and  the  transcription  end  site  (TES),  and 

309  not  within  the  promoter  range  of  the  same  gene.  To  assign  a  peak  to  an  intergenic  region  it  had  to  be  - 

310  10,000kb  from  the  TSS  and  -l-10,000kb  from  the  TES,  and  not  within  the  promoter  range  of  the  same  gene. 

311  A  peak  could  be  assigned  to  multiple  functional  regions  in  an  area  of  the  genome  harboring  multiple 

312  genes.  A  common  example  of  this  is  an  area  with  genes  on  both  strands.  A  peak  may  not  fit  any  of  these 

313  definitions  and  was  assigned  to  the  classification  “other”.  This  methodology  yielded  34,3 1 7  functional 

3 14  assignments  for  the  peaks  in  the  MC3T3-E1  cells. 

315  GEM  analysis:  Genome  wide  Event  finding  and  Motif  discovery  (GEM)  [Guo  et  ah,  2012]  was  used  to 

316  derive  the  Nmp4  consensus  sequence.  The  latest  mouse  genome  build  (mmlO)  was  employed  together 

317  with  the  GEM  default  ChlP-seq  read  distribution  file  and  a  minimal  k-mer  width  of  6  and  maximum  of  20. 

318  Gene  Ontology:  Gene  ontology  analysis  was  conducted  using  DAVID  [Huang  et  ah,  2009],  and  terms 

319  summarized  using  REVIGO  [Supek  et  ah,  2011].  The  ENCODE  ChlP-Seq  Significance  Tool  was 

320  employed  to  identify  enriched  transcription  factors  in  our  Nmp4  gene  target  list  [Auerbach  et  ah,  2013]. 

321  Additionally  some  functional  analysis  was  also  generated  through  the  use  of  QlAGEN’s  Ingenuity 

322  Pathway  Analysis  (IP A®,  QIAGEN  Redwood  City,  www.qiagen.com/ingenuity). 
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Bone  phenotype  statistical  analysis:  Statistical  evaluations  were  processed  using  the  program  JMP 
version  7.0.1  (SAS  Institute,  Cary,  NC).  The  animal  studies  employed  a  two-way  ANOVA  using 
genotype  and  treatment  as  the  independent  variables  followed  by  either  a  Tukey  HSD  or  LS  Means  post 
hoc  test  if  a  genotype  x  treatment  interaction  was  indicated.  Statistical  significance  was  set  at  p<0.05.  To 
compare  growth  rates  of  the  WT  and  Nmp4'^^  MSPCs  derived  from  various  experimental  mice,  we 
evaluated  the  slopes  of  log-transformed  cell  counts  regressed  onto  experimental  day  using  a  t-test.  The 
numbers  of  mice  per  treatment  group  and  replicates/treatment  for  the  cell  studies  are  indicated  in  the 
appropriate  figures  and  tables. 

RESULTS: 

Nmp4^'^  mice  are  not  protected  from  ovx-induced  bone  loss 

To  determine  whether  genetically  disabling  Nmp4  activity  protects  mice  from  ovx-induced  bone 
loss  as  it  does  from  unloading-associated  osteopenia  [Hino  et  al.,  2007],  we  removed  the  ovaries  or 
performed  sham  operations  on  both  WT  and  Nmp4  '^  mice  (Figure  1).  Both  the  ovx  WT  and  ovx  Nmp4^^~ 
mice  experienced  significant  weight  gain  at  4wks  post-op  (Table  1)  consistent  with  previous  mouse 
studies  [Vieira  Potter  et  al.,  2012].  Additionally,  ovx  resulted  in  a  significant  decrease  in  uterine  weight  in 
both  genotypes  (Table  1).  There  was  no  genotype  x  treatment  interaction  in  either  of  these  parameters. 

Both  WT  and  Nmp4^^^  mice  exhibited  significant  bone  loss  4wks  after  ovx  surgery  as  measured  in 
the  trabecular  bone  compartment  of  the  distal  femur  and  the  L5  vertebra  (Table  1).  The  Nmp4^'^  mice 
exhibited  a  trend  towards  enhanced  loss  of  bone  that  neared  significance  in  the  distal  femur  (BV/TV, 
genotype  x  treatment  interaction  =  0.06,  Table  1)  and  reached  significance  in  the  L5  vertebra  (BV/TV, 
genotype  x  treatment  interaction  <0.05,  Table  1).  Despite  this  enhanced  (or  nearly  enhanced)  rate  of  bone 
loss  the  Nmp4^^'  animals  maintained  more  trabecular  bone  compared  to  WT  mice  during  the  first  4wks 
after  ovariectomy.  The  Am/74  ^  mice  exhibited  a  decrease  in  the  serum  bone  formation  marker  PINP  at 
4wks  post-op  and  both  genotypes  showed  significant  decreases  in  this  marker  at  12wks  post-op  in  the 
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vehicle-treated  mice.  A  small  decrease  in  the  serum  bone  resorption  marker  CTX  was  observed  at  12wks 
post -op  in  the  vehicle-treated  mice. 

Ovx  Nmp4'^^  mice  show  an  enhanced  bone  gain  response  to  PTH  therapy 

With  a  separate  group  of  ovx  mice  we  initiated  treatment  of  both  WT  and  Nmp4-mA\  ovx  animals 
with  PTH  (30pg/kg/day)  or  vehicle  control  4wks  after  surgery.  The  duration  of  hormone  therapy  lasted 
4wks  (8wks  post-op)  and  8wks  (12wks  post-op).  The  ovx  Nmp4^'  mice  showed  an  enhanced  PTH- 
induced  gain  in  femoral  BV/TV  and  Conn  D  at  4wks  and  8wks  of  therapy  compared  to  their  ovx  WT 
littermates  as  well  as  an  augmented  gain  in  trabecular  thickness  at  8wks  (Figure  2,  Table  2).  The  null 
mice  also  showed  an  enhanced  PTH  response  at  the  L5  vertebra  at  8wks  of  treatment  (Figure  3,  Table  2). 
Specifically  the  2-way  ANOVA  indicated  strong  genotype  x  treatment  effects  for  the  distal  femur  for  both 
4wks  and  8wks  therapy  and  for  the  L5  vertebra  for  8wks  therapy  (see  Figures  2A  and  3 A);  the  post-hoc 
tests  concluded  that  the  difference  between  the  genotypes  was  within  the  hormone -treated  groups.  The 
vehicle-treated  ovx  WT  and  ovx  NmpP^'  groups  showed  no  difference  in  BV/TV  (Figures  2  and  3)  at  the 
end  of  the  treatment  regimens  indicating  that  the  modest  enhanced  loss  in  bone  in  the  Nmp4^^^  mice  was 
stabilized  by  4wks  therapy.  PTH  significantly  elevated  MAR,  MS/BS,  and  BFR  at  the  end  of  4wks 
treatment  as  shown  by  strong  treatment  effects  (Table  3).  However,  there  was  no  genotype  effect  or 
genotype  x  treatment  interaction  for  any  of  these  parameters  (Table  3).  Hormone  significantly  elevated 
serum  levels  of  the  bone  formation  marker  PI  NP  and  the  resorption  marker  CTX  at  8wks  of  therapy,  but 
there  was  no  treatment  x  genotype  interaction  for  either  of  these  parameters  (Table  3). 

FACS  analysis  of  the  BM  CD457CD105VCD146VnestiW  osteoprogenitors  revealed  a  significant 
elevation  in  the  number  of  these  cells  in  the  BM  obtained  from  the  Nmp4^^^  mice  at  the  end  of  4wk  therapy, 
irrespective  of  treatment  (Figure  4A).  This  is  consistent  with  our  previous  observation  in  the  ovary-intact 
null  mice  [He  et  ah,  2013].  By  the  end  of  8wks  treatment  (12wks  post-op)  the  observed  increase  in  the 
number  of  these  NmpP^'  cells  in  the  BM  failed  to  reach  statistical  significance,  but  there  was  a  significant 
elevation  in  the  number  of  the  PBL  NmpP^'  osteoprogenitors  in  the  vehicle-treated  mice  (Figure  4D).  The 
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375  Nmp4  '^  mice  showed  a  significant  elevation  in  CD8^  T  cells  in  both  the  BM  and  the  PBL  throughout  the 

376  entire  therapy  regimen  (Figure  4B  and  4E).  PTFl  significantly  decreased  the  numbers  of  these  cells  in  the 

377  BM  at  8wks  therapy  in  both  genotypes  (Figure  4B)  but  had  no  impact  on  the  number  of  these  cells  in  the 

378  PBL  (Figure  4E).  Disabling  Nmp4  had  little  to  no  effect  on  CD4^  T  cells,  nor  did  treatment  with  PTFl 

379  (Figure  4C  and  4F).  The  modest  increase  in  BM  CD4^  T  cells  approached  significance  (p<0.06)  but  this 

380  was  not  reflected  in  the  PBL,  just  as  we  previously  observed  in  the  ovary-intact  mice  [Fie  et  ah,  2013]. 

381  To  determine  if  the  enhanced  osteogenic  potential  of  the  BM  could  be  reliably  and  reproducibly 

382  maintained  in  vitro  in  MSPC  cultures  over  several  passages  and  in  the  absence  of  supporting  cells  (e.g.  T- 

383  cells)  we  established  expanded  WT  and  Nmp4  ^^  MSPCs  from  ovary-intact  mice.  The  expanded  Nmp4^^^ 

384  MSPCs  from  ovary-intact  mice  exhibited  modest  but  signiflcantly  enhanced  proliferation  compared  to  the 

385  WT  cells  (Figure  5A).  Both  the  null  and  WT  expanded  MSPCs  showed  strong  alkaline  phosphatase 

386  expression  (Figure  5B).  However,  the  expanded  Nmp4  '^  MSPCs  typically  showed  an  accelerated  and 

387  enhanced  mineralization  compared  to  WT  cells  under  various  concentrations  of  dexamethasone  and 

388  ascorbic  acid  (Figure  5B).  Finally,  the  expanded  Nmp4^^^  and  WT  MSPCs  exhibited  varying  degrees  of 

389  alkaline  phosphatase  staining  while  maintained  in  MesenCult  medium,  depending  on  the  confluence  of 

390  the  cells  and  time  in  culture  (3-9  days),  however  no  mineralization  was  observed  in  these  control  cultures 

391  (data  not  shown). 

392 

393  Genome-wide  ChIP-seq/gene  ontology  analysis  reveals  Nmp4  target  genes  and  potential  pathways  of  the 

394  anti-anabolic  axis. 

395  Nmp4  is  expressed  in  nearly  all  cells,  yet  the  most  singular  consequence  of  globally  disabling  this 

396  protein  is  the  enhanced  mobilization  of  bone  cells  upon  osteoanabolic  induction  [He  et  ah,  2013; 

397  Childress  et  ah,  2011;  Robling  et  ah,  2009;  Morinobu  et  ah,  2005].  As  a  first  step  in  understanding  the 

398  origins  of  this  phenotype,  which  may  have  clinical  significance,  we  needed  the  following  information:  (1) 

399  the  identity  of  the  Nmp4  target  genes  including  ‘core’  target  genes  common  to  multiple  cell  types;  (2) 

400  identify  common  functions  of  these  core  genes  to  distinguish  pathways  that  make  osteoprogenitors 
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401  particularly  vulnerable  to  the  effects  of  Nmp4  and  (3)  experimental  confirmation  of  some  of  these 

402  pathways.  To  begin  to  understand  how  Nmp4  works  we  set  out  to  understand  (4)  whether  Nmp4  targets 

403  functional  regions  of  the  genome,  (5)  if  it  binds  directly  to  DNA  or  via  other  proteins,  and  (6)  whether 

404  osteoanabolic  agents,  e.g.  PTH,  alter  Nmp4  DNA-binding  along  target  genes. 

405  The  potential  Nmp4  target  genes  identified  by  ChlP-seq  in  the  MC3T3-E1  (vehicle -treated)  cells 

406  and  those  established  in  the  three  ENCODE  cell  lines  were  compared  using  those  genes  that  had  one  or 

407  more  peaks  associated  with  the  TSS.  A  Venn  diagram  of  these  genes  showed  that  21 14  Nmp4  ‘core’ 

408  target  genes  were  common  to  the  four  cell  lines  (Figure  6A,  and  Supplemental  Table  1).  These  core  target 

409  genes  were  classified  into  functionally  related  categories  using  gene  ontology  (GO)  analysis  with  the 

410  Database  for  Annotation,  Visualization,  and  Integrated  Discovery  (DAVID)  tool  [Huang  et  al,  2009].  The 

411  functional  annotation-clustering  algorithm  was  applied  to  the  target  list,  which  is  able  to  give  a  more 

412  insightful  view  of  the  relationships  between  annotation  categories  and  terms  compared  to  other  analytic 

413  modules  [Huang  et  al.,  2009].  The  significance  of  group  classification  was  defined  by  enrichment  scores 

414  based  on  Fisher  exact  statistics  (false  discovery  rate,  FDR  p<0.05).  The  DAVID-derived  biological 

415  profile  was  further  summarized  using  REVIGO  [Supek  et  al.,  2011].  GO  analysis  of  the  core  target  genes 

416  designated  Nmp4  as  a  negative  regulator  of  cellular  biosynthetic  processes  showing  significant 

417  enrichment  for  genes  involved  in  the  regulation  of  transcription,  chromatin  modification,  protein  catabolic 

418  processes,  regulation  of  the  cell  cycle,  and  mRNA  processing/splicing  (Figure  6B).  Interestingly,  the 

419  genes  specific  to  any  one  particular  cell  line  or  specific  to  vehicle-treated  or  PTH-treated  MC3T3-E1  cells 

420  did  not  yield  a  distinct  biological  process  profile  that  reached  statistical  significance  as  obtained  with  the 

421  core  target  genes  (data  not  shown).  However,  peak-associated  genes  common  to  the  vehicle-  and  PTH- 

422  treated  MC3T3-E1  cells  yielded  a  profile  nearly  identical  to  that  obtained  with  the  core  target  genes. 

423  Next  we  probed  existing  datasets  for  enriched  transcription  factors  within  our  Nmp4  core  target 

424  gene  list  using  the  ENCODE  ChIP-seq  Significance  Tool  [Auerbach  et  al.,  2013]  (Table  4).  This  profile 

425  shows  that  Nmp4  binding  in  the  promoter  regions  of  its  target  genes  predominantly  co-occurs  with 

426  proteins  that  regulate  chromatin  organization  and  with  proteins  that  contribute  to  maintaining 
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stem/progenitor  pluripotency/multipotency  and  the  poised  gene  state,  e.g.  CHD2,  SIN3a,  and  GCN5 
[Harada  et  al.,  2012;  Nascimento  et  al.,  2011;  Lin  et  al.,  2007]. 

In  an  effort  to  gain  further  understanding  of  how  Nmp4  regulates  gene  expression  we  prepared  a 
genome -wide  functional  region  map  of  the  Nmp4  binding  sites  for  all  four  cell  types  as  described  in 
Materials  and  Methods.  The  majority  of  the  occupancy  peaks  were  located  in  or  near  the  TSS  or  in 
intragenic  regions,  areas  typically  associated  with  regulatory  functions  (Figure  7A).  To  determine  if 
Nmp4  binds  directly  to  DNA  or  can  associate  with  the  genome  via  other  proteins  we  used  the  discovery 
algorithm  GEM  to  derive  the  Nmp4  consensus-binding  site  from  the  MC3T3-E1  data.  In  support  of 
previous  studies  by  our  lab  and  others  the  derived  binding  site  matched  the  unusual  homopolymeric 
(dA-dT)  consensus  sequence  previously  derived  by  cyclic  amplification  and  electrophoretic  mobility  shift 
assay  [Alvarez  et  al.,  1998;  Nakamoto  et  al.,  2000]  (Figure  7B).  No  other  consensus  sequences  were 
identified  suggesting  a  single  and  direct  mode  of  genome  association,  mediated  by  the  Cys2Flis2  DNA- 
binding  domain  [Torrungruang  et  al.,  2002].  To  determine  whether  PTH  challenge  altered  Nmp4  DNA- 
binding  along  target  genes  we  generated  genome-wide  Nmp4  ChIP-seq  profiles  using  the  pre -osteoblast 
cell  line  MC3T3-E1  treated  with  hPTH(l-34)  or  vehicle  control  for  Ihr.  We  used  the  Ihr  time  point 
because  we  observed  the  most  significant  differences  in  femoral  mRNA  expression  profiles  between  WT 
and  Nmp4''^  mice  Ihr  after  injection  [Childress  et  al.,  2011].  Hormone  reduced  Nmp4  genome -wide 
occupancy  from  a  total  of  15,446  to  13,109  binding  sites.  However,  at  the  level  of  the  single  gene  there 
was  a  diversity  of  changes  in  Nmp4  occupancy,  i.e.  PTH  was  observed  to  remove  (e.g.  Nid2),  induce  (e.g. 
Ccdc53)  or  have  no  effect  on  Nmp4-DNA  association  (e.g.  Akt2,  Arrb2)  (Figure  8;  also  see  ChIP-qPCR 
confirmation  of  Nmp4  binding.  Figure  SI). 

As  a  first  step  in  the  validation  of  the  ChlP-seq-derived  anti-anabolic  map  we  interrogated  90 
mRNA  transcripts  in  non-differentiating  and  osteogenic  differentiating  WT  and  Nmp4'^'  expanded  MSPCs 
at  five  different  time  points.  The  accelerated  and  enhanced  mineralization  of  the  Nmp4^^^  MSPCs  (Figure 
5)  is  consistent  with  our  previous  observation  that  in  response  to  PTH  Nmp4^^^  mice  add  more  bone  and 
add  it  faster  than  WT  mice  [Childress  et  al.,  2011].  Our  choice  of  Nmp4  target  genes  (Supplemental  Table 
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453  1)  was  based  on  our  DAVID  analysis.  We  chose  both  core  target  genes  and  Nmp4  target  genes  identified 

454  in  the  MC3T3-E1  pre -osteoblasts.  DAVID  also  uses  KEGG  (Kyoto  Encyclopedia  of  Genes  and 

455  Genomes)  database  to  map  large  gene  lists  to  signaling  pathways  [Huang  et  ah,  2009].  For  example  the 

456  DAVID/KEGG  profile  of  the  Nmp4  core  target  genes  included  the  TOR  and  insulin/IGF  1  signaling 

457  pathways  (Table  5)  and  indeed  the  insulin/IGF l->IRSl->PI3K->PDKl->Akt  signaling  response  limb  is 

458  common  to  many  of  the  pathways  listed.  This  is  also  consistent  with  our  IPA  analysis  (Supplemental 

459  Table  3).  Also  included  were  Nmp4  target  genes  coding  for  proteins  involved  in  the  ubiquitin-proteasome 

460  system,  chromatin  remodeling,  transcription  regulation,  and  RNA  processing.  Finally  we  analyzed  the 

461  expression  of  osteogenic  differentiation  markers. 

462  Volcano  plots  (Supplemental  Figure  2)  identify  genes  that  were  statistically  significantly 

463  upregulated  or  downregulated  by  at  least  2-fold  in  the  Nmp4^^'  cells  compared  to  the  WT  cohort  on  the 

464  same  day  of  culture.  Figures  9  and  10  compare  the  relative  mRNA  expression  of  select  genes  to  the  level 

465  of  the  transcript  in  WT  cells  on  Day  3  of  culture  thus  providing  a  time  course  view.  The  Nmp4'^'  cells 

466  showed  a  strikingly  elevated  expression  of  the  Nmp4  target  gene  Igfbp2  mRNA  and  downregulation  of 

467  the  target  gene  Igfbp4  mRNA  during  the  early  differentiation  period  (Figures  9A  and  9B).  Igfbp2 

468  stimulates  osteoblast  differentiation  whereas  the  Igfbp4  is  a  potent  inhibitor  of  Igf  actions  [Xi  et  ah,  2014; 

469  Miyakoshi  et  ah,  1999].  Phosphoinositide-dependent  kinase  1  (Pdkl,  target  gene)  a  key  component  of  the 

470  Igfl/insulin  signaling  pathway  [Calleja  et  ah,  2014]  was  upregulated  in  the  NmpP^'  cells  throughout  the 

471  developmental  time  course  (Figure  9C).  Interestingly,  neither  the  target  genes  Igfl  nor  its  receptor  Igfl r 

472  exhibited  striking  differences  in  gene  expression  between  the  two  genotypes  (data  not  shown).  The  Nmp4' 

473  cells  exhibited  an  enhanced  anabolic  profde  during  the  latter  differentiation  period  as  evidenced  by 

474  elevated  expression  levels  of  the  non-target  genes  Bmp2  (Figure  9D),  Pthlr  (Figure  9E),  and  Bglap 

475  (osteocalcin.  Figure  9F). 

476  Cxcll2  expression  (target  gene),  also  known  as  stromal  derived  factor  1,  was  dramatically 

477  downregulated  in  the  NmpP^^  cells  throughout  development  (Figure  lOA)  and  the  target  gene  Plaur 

478  (uPAR,  urokinase  plasminogen  activator  receptor)  was  upregulated  in  the  null  cells  (Figure  lOB).  Both 
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genes  play  roles  in  MSPC  osteogenic  lineage  commitment  (Shahnazari  et  al.,  2013;  Sugiyama  et  al., 

2006;  Furlan  et  al.,  2007;  Kalbasi  et  al.,  2014).  Sppl  (osteopontin,  target  gene)  and  Thbs2 
(thrombospondin  2,  target  gene)  regulate  aspects  of  mineralization  [Hunter  2013;  Alford  et  al.,  2010]  and 
the  former  was  upregulated  in  our  null  cell  whereas  the  latter  was  downregulated  (Figures  IOC  and  lOD). 
Type  1  collagen  {Collal,  target  gene)  expression  was  elevated  in  the  Nmp4^^'  cells  throughout  the 
developmental  period  (Figure  lOE).  Interestingly,  we  observed  no  substantial  difference  in  the  expression 
profdes  of  Sp7  (osterix.  Figure  lOF)  or  the  target  gene  Runx2  (data  not  shown),  essential  transcription 
factors  for  osteoblast  differentiation  (Komori,  2011). 

DISCUSSION: 

Our  findings  that  the  Am/74  ^  mice  are  not  protected  from  ovarectomy-induced  bone  loss  yet 
maintain  the  amplified  response  to  PTH  therapy  is  a  key  advance  necessary  for  further  consideration  of 
Nmp4-based  treatment  strategies.  The  ovx  Nmp4'^  mice  displayed  an  enhanced  hormone -induced 
recovery  of  femoral  and  L5  trabecular  BV/TV  despite  delaying  treatment  until  4wks  post-op  to  allow  for 
significant  bone  loss.  Both  the  ovx  WT  and  ovx  Nmp4  '^  mice  showed  strong  responses  to  PTH  therapy. 
After  4wks  and  8wks  of  treatment  the  WT  mice  displayed  a  3.2-fold  and  4.6-fold  increase  in  femoral 
BV/TV  over  vehicle -treated  mice,  respectively.  However  the  Nmp4^'^  mice  showed  a  3.6-fold  and  8.8-fold 
increase  over  the  same  time  period  resulting  in  a  very  strong  genotype  x  treatment  interaction.  Differences 
in  PTH-mediated  BV/TV  restoration  efficacy  between  the  WT  and  Nmp4^^^  mice  in  the  L5  vertebra  and 
was  less  striking  although  statistically  significant  (1.3 -fold  vs  1.6-fold  at  8wks  in  the  WT  and  Nmp4'^' 
mice,  respectively).  We  observed  similar  PTH-responsive  femoral  and  L5  profiles  between  younger, 
ovary-intact  WT  and  Nmp4-mi\\  mice  [Robling  et  al.,  2009;  Childress  et  al.,  201 1;  He  et  al.,  2013].  The 
histomorphometry  and  serum  data  reported  here  tracked  the  PTH-induced  increases  in  bone  mass  in  the 
ovx  animals  showing  strong  treatment  effects  for  bone  formation  parameters  MAR,  BFR,  and  MS/BS  (at 
4wks  treatment)  as  well  as  strong  increases  in  bone  remodeling  serum  PI NP  and  CTX  (at  8wks  treatment). 
However,  these  parameters  did  not  distinguish  the  genotypes  in  regards  to  the  amount  of  bone  formed 
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over  this  time  period  as  was  achieved  with  the  pCT  data.  Interestingly,  the  histomorphometry  data  did  not 
distinguish  the  differences  in  PTH-induced  bone  formation  in  ovary-intact  WT  and  Nmp4^^^  mice 
[Childress  et  ah,  201 1],  Our  present  observation  that  Nmp4^^'  MSPCs  exhibit  an  accelerated  and  enhanced 
mineralization  in  response  to  anabolic  cues,  i.e.  osteogenic  medium  suggests  that  the  augmented  bone 
formation  is  an  early  event.  Similarly,  we  did  not  observe  the  expected  ovx-induced  small  increase  in 
serum  CTX.  Instead,  the  serum  data  of  the  vehicle-treated  mice  showed  a  large  decrease  in  PINP  and  a 
smaller  decrease  in  CTX  over  the  time  course  of  the  entire  experiment,  i.e.  pre-op  vs  12wks  post-op.  This 
suggests  a  decrease  in  bone  remodeling  in  the  untreated  mice  predominantly  in  the  bone  formation  arm.  A 
more  extensive  time  course  with  earlier  harvest  points  for  histomorphometry  and  serum  samples  is 
required  to  more  fully  characterize  the  anticipated  differences  in  WT  and  Nmp4'^’  dynamic  bone 
remodeling. 

The  most  robust  phenotypic  characteristic  of  Nmp4  ablation  is  the  exaggerated  bone  formation 
response  to  PTH  or  BMP2,  which  suggests  that  the  adult  mice  harbor  an  increased  number  of  BM  MSPCs 
with  heightened  sensitivity  to  osteoanabolic  signals.  Disabling  Nmp4  has  no  observable  impact  on 
embryonic  or  perinatal  skeletal  development.  Adult  MSPCs  are  a  heterogeneous  population  of 
multipotent  stem,  progenitor,  and  stromal  cells  that  contribute  to  BM  homeostasis  [Mizoguchi  et  ah, 
2014].  In  mouse  bone  marrow  much  of  the  CFU-F  activity  is  in  the  nestiW  cell  population  and  in  the 
human  marrow  the  CD146^  population  [Mizoguchi  et  ah,  2014;  Sacchetti  et  ah,  2007].  In  ovary-intact, 
Nmp4  '^  mice  we  observed  a  4-fold  increase  in  the  frequency  of  CD45VCD105TCD146VnestiW  cells 
irrespective  of  treatment  (PTH  vs  vehicle  control),  which  paralleled  the  magnitude  increase  in  CFU-F  and 
CFU-F^"^P'*°'‘^  cell  number  in  culture  [He  et  ah,  2013].  Similarly,  the  ovx  Nmp4~^'  mice  exhibited  an 
approximate  3 -fold  increase  in  the  CD457CD105TCD146AnestiW  cells  at  8wks  post-op  compared  to  the 
ovx  WT  animals. 

The  enhanced  osteogenic  potential  of  the  BM  as  measured  by  the  frequency  of  cells 

capable  of  becoming  osteoprogenitors  persists  in  expanded  Nmp4^^^  MSPC  cultures  over  5-10  passages 
and  removed  from  the  supporting  CD8^  T  cells.  In  culture  these  cells  displayed  a  modest  increase  in 
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proliferative  activity  and  perhaps  this  aspect  of  the  phenotype  contributes  to  the  observed  expanded  pool 
of  osteoprogenitors  in  vivo.  In  an  earlier  study,  Noda  and  colleagues  demonstrated  that  Nmp4^^  BM 
yielded  significantly  more  CFU-F®’’  mineralizing  colonies  at  passage  Po  than  WT  BM  [Morinobu  et  ah, 
2005].  Our  present  data  extend  these  observations  and  show  that  the  serially  passaged  Nmp4-/-  MSPCs 
maintain  a  strikingly  enhanced  capacity  for  mineralization  compared  to  the  capacity  of  the  WT  cultures. 
Taken  together  these  observations  suggest  that  there  is  a  cell  autonomous  role  of  Nmp4  for  regulating 
MSPC  osteogenesis. 

Our  -omics  data  combined  with  our  low-density  array  results  suggest  that  upon  challenge  with  an 
anabolic  cue  Nmp4'^^  MSPCs  produce  autocrine/paracrine  factors  that  enhance  the  replication  and 
differentiation  of  neighboring  osteoprogenitors,  a  key  early  event  driving  the  PTH  anabolic  response 
[Jilka,  2007].  We  observed  that  once  Nmp4'^'  cells  were  transferred  to  osteogenic  medium  they  expressed 
strikingly  elevated  levels  of  the  Nmp4  target  gene  Igfbpl,  a  strong  autocrine/paracrine  factor  that 
enhances  osteogenesis  (Xi  et  ah,  2014).  Consistent  with  our  observations,  overexpression  of  Igfbp2  in 
MC3T3-E1  cells  accelerated  the  time  course  of  differentiation  and  mineralization  as  well  as  increased  the 
total  number  of  differentiating  cells.  By  Day  6  in  this  previous  study,  Igfbp2-overexpressing  cells 
expressed  twice  as  much  osteocalcin  as  control  cultures  and  this  difference  persisted  [Xi  et  ah,  2014]. 

This  is  a  strikingly  similar  phenotype  to  the  Nmp4'^'  cells.  Interestingly,  the  expression  of  the  Nmp4  target 
gene  lgfbp4  was  decreased  in  the  null  cells.  This  binding  protein  is  a  strong  inhibitor  of  osteoblast 
differentiation  [Miyakoshi  et  ah,  1999]  and  thus  its  suppression  may  further  accelerate  and  enhance  the 
differentiation  of  the  null  cells.  Igfl  is  a  key  mediator  of  the  PTH  anabolic  response  [Bikle  &  Wang, 

2012;  Elis  et  ah,  2010]  and  although  there  was  no  notable  alteration  in  the  expression  profiles  of  Igfl  or 
Igflr  'vn  the  null  cells,  the  expression  of  Pdkl,  a  target  of  Nmp4  and  a  key  kinase  component  of  the 
Igfl/insulin  signaling  pathway  was  elevated.  This  may  enhance  the  sensitivity  of  the  Nmp4'^'  cells  to  this 
growth  factor.  At  the  end  of  the  16-day  culture  period  the  Nmp4'^'  cells  exhibited  an  enhanced  anabolic 
profile  as  evidenced  by  the  elevated  expressions  of  the  non-target  Nmp4  genes  Bmp2,  Plhlr,  and  Bglap. 
Perhaps  this  is  an  autocrine/paracrine  response  to  the  earlier  surge  in  Igfbpl  expression.  Indeed  the  Igfl 
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557  pathway  plays  a  significant  role  in  MSPC  proliferation  and  mineralization  [Kumar  and  Ponnazhagan, 

558  2012;  Xian  et  al.,  2012]  and  the  null  cells  exhibited  alterations  in  the  expression  not  only  of  Bglap  but  the 

559  Nmp4  target  extracelluar  matrix  proteins  Collal,  Sppl,  and  Thbsl.  Although  the  molecular  mechanisms 

560  underlying  mineralization  remain  to  be  elucidated,  Sppl  is  an  anionic  phosphoprotein  expressed  in 

561  mineralizing  tissues  that  appears  to  regulate  crystal  size,  shape,  and  location  [Hunter,  2013].  Thbsl  is  an 

562  extracellular  matrix  glycoprotein  that  has  pleiotropic  effects  on  bone  phenotype.  This  protein  appears  to 

563  suppress  the  MSPC  osteoprogenitor  pool  but  also  supports  mineralization  [Alford  et  al.,  2010;  Hankenson 

564  et  al.,  2002;  Hankenson  et  al.,  2000].  Therefore  whether  the  decrease  in  Thbsl  expression  in  the  NmpT^^ 

565  cells  impacts  the  observed  alteration  in  the  number  of  osteoprogenitors,  alterations  in  mineralization  or 

566  impacts  cell  phenotype  in  other  ways  remains  to  be  determined.  Finally,  we  observed  no  striking 

567  differences  in  the  expression  of  the  transcription  regulators  Sp7  and  Runxl  between  the  genotypes.  This 

568  suggests  that  disabling  Nmp4  alters  select  aspects  of  the  developing  osteoblast  phenotype. 

569  The  dramatic  decrease  in  Cxclll  expression  in  the  Nmp4'^'  cells  raises  the  question  as  to  whether 

570  this  plays  a  role  in  the  observed  increase  in  osteoprogenitors  in  the  null  mice  [He  et  al., 

571  2013;  Morinobu  et  al.,  2005  and  the  present  work].  CxcH2  and  its  receptor  Cxcr4  play  key  roles  in 

572  maintaining  the  bone  marrow  niche  and  CxcH2  is  expressed  by  bone  marrow  stromal  cells  and  cells  of 

573  the  osteoblast  lineage  [Jung  et  al.,  2006;  Sugiyama  et  al.,  2006].  Ablation  of  the  receptor  Cxcr4  in  mature 

574  osteoblasts  increased  the  number  of  cFU-F'‘"‘p*'°'®^  osteoprogenitors  recovered  from  these  mice  although 

575  the  phenotype  also  included  a  decrease  in  BV/TV  [Shahnazari  et  al.,  2013].  Our  results  suggest  that 

576  suppressing  the  expression  of  the  Cxcr4  ligand  CxcH2  results  in  a  similar  impact  on  osteoprogenitor 

577  number  but  a  different  bone  phenotype.  The  upregulation  of  Plaur  expression  in  the  Nmp4~^^  MSPCs  may 

578  potentially  contribute  to  the  increased  number  of  CFU-F‘‘'''p''°^^  cells  since  abrogating  the  activity  of  this 

579  GPI-anchored  receptor  suppressed  MSPC  osteogenic  differentiation  [Kalbasi  et  al.,  2014].  Finally 

580  mutagenesis  and  rescuing  experiments  to  determine  whether  the  potential  targets  are  truly  functionally 

581  significant  is  the  next  required  step  for  authenticating  this  new  anti -anabolic  network. 
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Further  parsing  of  the  enhanced  Nmp4''^  BM  osteogenic  potential  implicates  the  elevated 
frequency  of  CD8^  T  cells  in  both  ovary-intact  and  ovx  Nmp4^^^  mice,  although  this  requires  functional 
confirmation  in  these  models.  The  ovx  null  animals  exhibited  elevated  numbers  of  CD8^  T  cells  in  both 
BM  and  PBL  compartments  throughout  the  entire  treatment  regimen,  similar  to  what  we  previously 
observed  in  the  younger  ovary-intact  Nmp4  '^  mice,  although  this  increase  was  limited  to  the  BM  [He  et  ah, 
2013].  The  elevated  number  of  CD8^  T  cells  is  intriguing  since  these  cells  are  documented  to  amplify  the 
PTH  anabolic  response  [Bedi  et  al.,  2012;  Terauchi  et  al.,  2009].  MSPCs  regulate  T  cell  proliferation  and 
survival  [Wang  et  al.,  2012]  and  perhaps  disabling  Nmp4  de -represses  this  aspect  of  the  cell-cell 
interaction,  although  this  apparent  alteration  in  proliferation/survival  may  be  a  cell  autonomous  feature  of 
the  Nmp4'^'  T  cell  phenotype.  Although  the  elevated  number  of  Am/74'^  osteoprogenitors  declined  at  12wks 
post -op,  the  higher  number  of  Nmp4'^^  CD8^  T  cells  did  not  decrease.  Therefore,  although  the  augmented 
response  to  PTH  might  be  weakened  in  the  Nmp4‘^'  mice,  it  may  not  disappear  since  the  persistent 
increased  lymphocyte  number  might  provide  extra  WntlOb  as  a  potent  osteoprogenitor  differentiation 
factor. 

Why  is  there  typically  no  difference  between  the  amount  of  baseline  trabecular  bone  in  WT  and 
Nmp4  ''  mice  despite  the  presence  of  an  expanded  pool  of  osteoprogenitors  and  CD8^  T  cells  in  the  null 
bone  marrow?  This  phenomenon  differs  from  the  results  of  recent  clinical  trials  in  which  neutralizing 
sclerostin,  an  inhibitor  of  the  Wnt  signaling  pathway  and  osteoblast  differentiation,  significantly  increases 
baseline  bone  mineral  density  [reviewed  in  Becker,  2014].  Apparently  disabling  Nmp4  is  not  sufficient 
for  driving  excess  bone  formation  but  instead  primes  the  aforementioned  cells  for  activation  by  an 
anabolic  cue.  The  occasionally  observed  elevated  trabecular  volume  in  untreated  Nmp4^^^  mice  may  be 
due  to  the  sporadic  local  release  of  growth  factors,  e.g.,  Igfl,  lgfbp2  or  Bmp2.  However  exogenous 
pharmacological  doses  of  PTH  provide  a  strong  stimulus  for  triggering  the  response  leading  to  the 
enhanced  bone  formation.  Once  activated  by  the  anabolic  cue,  the  Nmp4^^'  cells  produce  the 
autocrine/paracrine  factors  that  enhance  the  anabolic  response. 
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607  Consistent  with  the  requirement  for  a  strong  anabolic  cue  to  trigger  enhanced  bone  formation  in 

608  the  Nmp4  '^  mice,  disabling  this  transcription  factor  did  not  protect  the  animals  from  ovx -induced  bone 

609  loss,  indeed  the  initial  rate  of  loss  during  the  first  4wks  after  ovariectomy  was  higher  (L5)  or  nearly 

610  higher  (distal  femur)  in  the  Nmp4^^^  mice.  These  animals  harbor  a  modestly  elevated  number  of  osteoclast 

611  progenitors  (CFU-GM)  [He  et  ah,  2013]  that  upon  differentiation  exhibit  an  enhanced  bone -resorbing 

612  activity  in  vitro  [Childress  et  ah,  201 1].  Therefore  a  decrease  in  estrogen  might  accentuate  this  aspect  of 

613  the  phenotype.  Moreover,  differences  in  sex  steroid  levels  may  underlie  why  intact  male  Nmp4^^^  mice  did 

614  not  lose  bone  under  hind  limb  suspension  [Hino  et  ah,  2007].  As  mentioned,  the  ^  baseline 

615  phenotype  includes  an  occasional  unprovoked  enhancement  in  trabecular  architecture,  which  we  observed 

616  in  the  present  study.  That  is  to  say,  despite  the  elevated  initial  bone  loss,  the  cohort  of  sham  and  ovx 

617  Nmp4  '^  mice  had  more  femoral  and  L5  trabecular  bone  compared  to  WT  at  the  time  of  harvest  (Table  3). 

618  However,  there  was  no  statistical  difference  between  vehicle-treated  animals  in  either  the  4wk  or  8wk 

619  hormone  therapy  cohorts  (Figures  2  and  3).  Longitudinal  studies  for  serum  turnover  markers  coupled  with 

620  pQCT  in  live  mice  could  be  used  to  track  the  real-time  dynamics  of  ovx -induced  bone  loss  and 

621  subsequent  therapy-induced  bone  gain  between  the  WT  and  Nmp4  ^^  mice.  In  lieu  of  this,  we  employed  a 

622  2-way  ANOVA,  which  incorporates  differences  in  control  groups,  to  evaluate  whether  there  is  an 

623  interaction  between  genotype  and  treatment. 

624  The  present  data  also  contribute  to  our  knowledge  as  to  how  Nmp4  works  at  the  molecular  level. 

625  Nmp4  binds  throughout  the  genome  but  is  primarily  localized  to  regions  near  the  TSS  and  within  the  gene, 

626  consistent  with  mediating  a  regulatory  role.  GEM  analysis  confirmed  the  AT -rich  homopolymeric 

627  binding-site  and  did  not  identify  other  consensus  sequences  expected  only  if  Nmp4  also  interacted  with 

628  the  genome  indirectly  via  other  DNA-binding  proteins.  Nmp4  association  with  the  genome  is  responsive 

629  to  PTH  since  hormone  decreased  genome -wide  occupancy  in  the  MC3T3-E1  cells  after  Ihr  of  exposure. 

630  However,  the  impact  of  PTH  on  Nmp4  occupancy  was  gene  and  site -specific  and  hormone  stimulation 

631  was  observed  to  induce,  remove,  or  have  no  effect  on  Nmp4  genomic  occupancy.  This  may  further 
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632  augment  the  fine  control  that  this  transcription  factor  has  over  the  regulation  of  osteoprogenitor  and/or 

633  bone-forming  capacity. 

634  There  is  a  critical  need  for  osteoanabolic  agents  [Lewiecki,  2011].  We  have  taken  a  two-pronged 

635  approach  in  our  research  to  serve  this  clinical  demand:  (1)  identify  molecular  and  cellular  mechanisms 

636  that  could  be  used,  for  example  in  an  adjuvant  setting  to  promote  enhanced  efficacy  or  less  frequent 

637  dosing  with  current  osteoanabolic  agents;  and  (2)  identify  innovative  approaches  to  identify  new  drug 

638  targets/pathways  or  mechanisms  of  action  that  would  provide  needed  substrate  for  the  future  drug 

639  discovery  initiatives  in  bone  disease,  including  osteoporosis.  Our  discovery-driven  approaches  have 

640  mapped  a  global  network  of  Nmp4-regulated  pathways  potentially  comprising  a  bone  anti-anabolic  axis. 

641  Further  functional  studies  charting  the  hierarchy  and  interactions  of  theses  network  pathways  will  provide 

642  a  novel  integrated  mechanism  underlying  the  natural  constraints  on  bone  formation.  We  postulate  that  the 

643  Nmp4  anti-anabolic  network  may  constitute  a  novel  strategy  to  identify  and  reveal  pharmacologically 

644  accessible  pathways  for  adding  new  bone  to  the  old  skeleton. 

645 
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897  FIGURE  LEGENDS 

898 

899  Eigure  1:  Schematic  of  treatment  regimen  for  WT  and  Nmp4^^'  mice;  Group  1  mice  were  subjected  to 

900  ovariectomy  (ovx)  or  sham  operation  at  12wks  of  age  and  evaluated  for  bone  loss  4wks  post-op  (16wks  of 

901  age).  Group  2  mice  were  ovx  at  12wks  of  age  and  began  PTH  or  vehicle  therapy  at  16wks  of  age  for  a 

902  duration  of  4wks  and  8wks.  Endpoint  analyses  included  micro-computed  tomography  pCT,  serum 

903  analysis  for  N-terminal  propeptide  of  type  1  procollagen  (PINP)  and  C-terminal  telopeptides  (CTX),  and 

904  dynamic  histomorphometry. 

905 

906  Eigure  2:  Disabling  Nmp4  enhances  PTH  restorative  therapy  in  the  distal  femur  of  ovx  Nmp4^^^  mice  [A] 

907  Interaction  plots  of  femoral  trabecular  bone  volume/total  volume  (BV/TV)  of  ovx  WT  and  ovx  Nmp4'^^ 

908  mice  as  determined  by  pCT  at  4wks  of  treatment  and  8wks  of  treatment.  Data  are  average  ±  SD,  number 

909  of  mice/experimental  group  =  8-9).  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

910  significance  was  set  at  p<0.05.  The  Tukey’s  HSD  post  hoc  test  was  used  to  determine  differences  between 

911  the  treatment  groups.  There  were  genotype,  treatment  and  genotype  x  treatment  interaction  at  both  time 

912  points.  There  was  no  difference  between  the  vehicle-treated  WT  and  Nmp4^^'  mice.  [B]  pCT  images 

913  showing  PTH-induced  improvements  in  distal  femur  trabecular  architecture  in  ovx  WT  and  Nmp4  ^'  mice 

914  after  8  weeks  of  treatment  (12wks  post-op,  24wks  of  age). 

915 

916  Eigure  3:  The  exaggerated  response  to  anabolic  PTH  persists  in  the  L5  vertebra  of  ovx  Nmp4^^^  mice.  [A] 

917  Interaction  plots  of  L5  vertebra  bone  volume/total  volume  (BV/TV)  of  ovx  WT  and  ovx  Nmp4^^'  mice  as 

918  determined  by  pCT  at  4wks  of  treatment  and  8wks  of  treatment.  Data  are  average  ±  SD,  number  of 

919  mice/experimental  group  =  8-9).  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

920  significance  was  set  at  p<0.05.  The  LS  Means  Student  t  post  hoc  test  was  used  to  determine  differences 

921  between  the  treatment  groups.  There  were  genotype,  treatment  effects  at  both  time  points  and  a  genotype 

922  X  treatment  interaction  at  8wks  therapy.  There  was  no  difference  between  the  vehicle-treated  WT  and 
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923  Nmp4  '^  mice.  [B]  |iCT  images  showing  PTH-induced  improvements  in  L5  trabecular  architecture  in  ovx 

924  WT  and  Nmp4  '^  mice  after  8  weeks  of  treatment  (12wks  post-op,  24wks  of  age). 

925 

926  Figure  4:  Ovx  does  not  abrogate  the  expanded  population  of  osteoprogenitors  and  CD8+  T  cells  in  Nmp4 

927  ^  mice.  FACS  analysis  of  BM  and  PBL  osteoprogenitors,  CD8+  T  cells,  and  CD4-I-  T  cells.  [A,  D]  The 

928  frequency  of  femoral  BM  and  PBL  CD45-/CD105+/CD146+/CD105+/nestin-l-  osteoprogenitor  cells  in 

929  WT  and  Nmp4  '^  mice  at  the  end  of  4wks  and  8wks  treatment  with  intermittent  PTH  or  vehicle  control;  [B, 

930  E]  the  frequency  of  BM  and  PBL  CD8+  T  cells  from  the  WT  and  Amp4  ^  mice;  [C,  F]  the  frequency  of 

931  BM  and  PBL  CD4-I-  T  cells  from  the  WT  and  null  mice.  Data  are  average  ±  SD,  number  of 

932  mice/experimental  group  =  8-9;  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

933  significance  was  set  at  p<0.05. 

934 

935  Figure  5:  Expanded  Nmp4  '^  MSPCs  exhibit  enhanced  proliferation  and  mineralization  in  culture.  [A] 

936  Comparative  growth  rates  of  expanded  WT  and  Nmp4  '^  MSPCs.  Cell  counts/day  (n=4  lines  per  genotype 

937  loglO  cells/well,  3  wells/sample,  average  ±  SD,  t  test,  t<0.05).  Note:  each  Tine’  is  derived  from  a  single 

938  mouse  [B]  Alkaline  phosphatase  (alk  phos)  and  alizarin  red  staining  of  a  WT  and  Amp4  ^  MSPC  cultures 

939  from  Day7-Day28.  See  text  for  details 

940 

941  Figure  6:  Nmp4  associates  with  core  target  genes  common  to  multiple  cell  types  and  acts  as  a  negative 

942  regulation  of  cellular  biosynthetic  processes  [A]  Venn  diagram  illustrating  the  shared  Nmp4  target  genes 

943  in  the  MC3T3-E1  osteoblast-like  cells  (vehicle-treated),  and  the  three  ENCODE  cells  lines,  ES-E14 

944  (embryonic  stem  cells),  MEL,  and  CH12  cells  (B-cell  lymphomas).  [B]  DAVID/REVIGO  gene  ontology 

945  (GO)  profile  of  Nmp4  core  target  genes 

946 

947  Figure  7:  Nmp4  binds  to  AT-rich  DNA  typically  proximal  to  TSS  sites  or  within  intragenic  regions.  [A] 

948  Genome -wide  mapping  of  the  Nmp4  binding  sites  show  that  most  sites  are  distributed  in  the  TSS  and 
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intragenic  regions  of  the  genome.  ChIP-seq  analysis  included  vehicle-treated  and  PTH-treated  MC3T3-E1 
osteoblast-like  cells  (vMC  and  pMC,  respectively)  and  three  murine  cell  lines  from  the  ENCODE 
Consortium  including  ES-E14  (Esl4),  which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and 
two  mouse  erythroleukemia  cell  lines  (Chl2  and  MEE)  derived  from  B-cell  lymphomas.  [B]  GEM 
analysis  for  the  Nmp4  consensus  sequence  derived  from  MC3T3-E1  cells.  A  minimal  k-mer  width  of  6 
and  maximum  of  20  were  used.  The  optimal  position  weight  matrix  (PWM)  score  for  the  MC3T3-E1  data 
was  10.07.  The  hypergeometric  P-value  (hgp)  was  le-1466.1. 

Figure  8:  ChlP-seq  reveals  Nmp4  binding  profiles  at  specific  gene  loci.  Mouse  MC3T3-E1  cells  were 
seeded  into  twenty-one  150mm  plates  at  an  initial  density  of  50,000  cells/plate  (320  cells/cm^)  and 
maintained  in  aMEM  complete  medium  -t  ascorbic  acid  for  14  days.  Prior  to  harvest  cells  were  treated 
with  25nM  hPTE[(l-34)  or  vehicle  control  for  Ihr.  Processing  for  ChIP-seq  analysis  was  performed  as 
described  in  the  Materials  and  Methods.  Sequences  (50nt  reads,  single  end)  were  aligned  to  the  mouse 
genome  (mm  10)  using  the  BWA  algorithm.  Alignments  were  extended  in  silico  at  their  3 ’-ends  to  a 
length  of  150bp,  which  is  the  average  genomic  fragment  length  in  the  size-selected  library,  and  assigned 
to  32-nt  bins  along  the  genome.  Nmp4  (Znfi84)  peak  locations  were  determined  using  the  MACS 
algorithm  (v  1.4.2)  with  a  cutoff  of  pvalue  =  le-7.  The  genomic  loci  including  the  chromosome  number 
and  nucleotide  interval  are  indicated.  Read  scales  are  indicated  on  the  Y-axis.  An  arrow  indicates  the 
transcriptional  start  sites  and  direction  of  transcription  for  each  of  the  genes;  vertical  boxes  within  the 
gene  indicate  exons.  The  Nmp4  ChIP-seq  gene  profiles  include  (A)  Nidi  (B)  Aktl,  (C)  Pdkl  (D)  ccdc53, 
(E)  Arrbl  and  (F)  Irsl.  The  input  DNA  profiles  were  devoid  of  peaks. 

Figure  9:  Comparison  of  mRNA  expression  profiles  derived  from  non-differentiating  (Day  3)  and 
osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^'  cells.  All  transcript  levels  are  compared  to  WT 
Day  3  providing  a  time  course  of  expression.  mRNA  profiles  [A]  Igfbpl;  [B]  lgfbp4;  [C]  Pdkt,  [D] 

Bmpl;  [E]  Pthlr  were  derived  from  the  TEDA  system  (Format  96a,  Applied  Biosystems,  Foster  City, 
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CA)  performed  on  a  QuantStudio^'^  7  Flex  Real-Time  PCR  System  and  normalized  with  GusB.  Profde 
[F]  Bglap  mRNA  profde  qRT-PCR  reactions  were  performed  on  an  Eppendorf  Mastercycler®  RealPlex^ 
using  Rplp2  as  the  normalizer  as  previously  described  [Robling  et  al,  2009],  Comparison  of  profdes  using 
GusB  and  Rplp2  as  the  normalizer  showed  no  differences  in  the  shape  of  the  expression  profiles 

Figure  10:  Comparison  of  mRNA  expression  profdes  derived  from  non-differentiating  (Day  3)  and 
osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^'  cells.  All  transcript  levels  are  compared  to  WT 
Day  3  providing  a  time  course  of  expression.  mRNA  profiles  [A]  Cxcll2;  [B]  Plaur;  [C]  Sppl;  [D] 

Thbs2;  [E]  Collal  were  derived  from  the  TEDA  system  (Format  96a,  Applied  Biosystems,  Foster  City, 
CA)  performed  on  a  QuantStudio^"^  7  Flex  Real-Time  PCR  System  and  normalized  with  GusB.  Profde 
[F]  Sp7  mRNA  profde  qRT-PCR  reactions  were  performed  on  an  Eppendorf  Mastercycler®  RealPlex^ 
using  Rplp2  as  the  normalizer  as  previously  described  [Robling  et  al.,  2009].  The  Day  16  WT  sample  is 
the  average  of  two  replicates.  Comparison  of  profiles  using  GusB  and  Rplp2  as  the  normalizer  showed  no 
differences  in  the  shape  of  the  expression  profdes 

Supplemental  Figure  SI:  qRT-PCR  validates  the  ChIP-seq  profiles.  [A]  The  Nmp4  ChlP-seq  profile  for 
the  gene  Collal.  The  genomic  loci  including  the  chromosome  number  and  nucleotide  interval  are 
indicated.  Read  scale  is  indicated  on  the  Y-axis.  An  arrow  marks  the  transcriptional  start  site  and  direction 
of  transcription;  vertical  boxes  within  the  gene  identify  exons.  [B]  qRT-PCR  was  used  to  authenticate  the 
ChIP-seq  peaks  as  described  in  the  Materials  and  Methods. 

Supplemental  Figure  2:  Volcano  plots  derived  from  gene  expression  profdes  of  non-differentiating  (Day 
3)  and  osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^'  as  described  in  the  Materials  and  Methods. 
RNA  expression  profiling  was  performed  on  a  QuantStudio^"^  7  Flex  Real-Time  PCR  System  and  data 
analyzed  using  the  ExpressionSuite  vl.O.d'^'^  analysis  software  (Applied  Biosystems)  as  described  in  the 
Materials  and  Methods.  [A]  WT  vs  Nmp4'^'  cells  at  Day  3  post-seeding.  mRNA  transcript  expression  was 
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compared  to  WT  cells  (Day  3).  Cells  maintained  in  Mesencult™  Media  +  Mesencult™  Stimulatory 
Supplement.  [B]  WT  vs  Nmp4'^'  cells  at  Day  5  post-seeding.  mRNA  transcript  expression  was  compared 
to  WT  cells  (Day  5).  Cells  maintained  in  differentiation  medium  for  48hrs  [C]  WT  vs  Nmp4'^'  cells  at  Day 
7  post-seeding.  mRNA  transcript  expression  was  compared  to  WT  cells  (Day  7).  Cells  maintained  in 
differentiation  medium  for  96hrs  [D]  WT  vs  Nmp4'^'  cells  at  Day  9  post-seeding.  mRNA  transcript 
expression  was  compared  to  WT  cells  (Day  9).  Cells  maintained  in  differentiation  medium  for  144hrs  [E] 
WT  vs  Nmp4''^  cells  at  Day  1 6  post-seeding.  mRNA  transcript  expression  was  compared  to  WT  cells 
(Day  16).  Cells  maintained  in  differentiation  medium  for  192hrs.  Genes  indicated  with  green  dots  (left  of 
center)  and  above  the  X-axis  exhibited  a  significant  downregulation  by  over  2-fold.  Genes  indicated  with 
the  red  dots  (right  of  center)  and  above  the  X-axis  exhibited  a  significant  upregulation  by  over  2-fold. 

Supplemental  TABLE  1  :  96  Nmp4  ‘core’  target  genes,  non-core  target  genes,  and  non-target  genes 
including  5  candidate  normalizer  genes.  Individual  cDNAs  were  quantified  by  qRT-PCR  using  a  custom 
TLDA  system  (Format  96a,  Applied  Biosystems,  Foster  City,  CA)  as  described  in  the  Materials  and 
Methods  on  a  QuantStudio™  7  Flex  Real-Time  PCR  System.  We  used  the  ExpressionSuite  vl.0.4™ 
analysis  software  (Applied  Biosystems)  to  analyze  these  data. 

Supplemental  TABLE  2:  21 14  Nmp4  ‘core’  target  genes  common  to  the  four  cell  lines  MC3T3-E1 
osteoblast-like  cells,  and  three  murine  cell  lines  from  the  ENCODE  Consortium  including  ES-E14  (Esl4), 
which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and  two  mouse  erythroleukemia  cell  lines 
(Chl2  and  MEE)  derived  from  B-cell  lymphomas. 

Supplemental  TABLE  3:  IPA  analysis  of  21 14  Nmp4  ‘core’  target  genes  common  to  the  four  cell  lines 
MC3T3-E1  osteoblast-like  cells,  and  three  murine  cell  lines  from  the  ENCODE  Consortium  including 
ES-E14  (Esl4),  which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and  two  mouse 
erythroleukemia  cell  lines  (Chl2  and  MEE)  derived  from  B-cell  lymphomas. 
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Table  1:  Bone  loss  data. 


WT 

Nmp4'^' 

2-WAY  ANOVA  n-values 

SHAM  OVX 

SHAM  OVX 

Genotype 

Treatment 

Gene  x  Treat 

%A  Body  weight 

2.48±7.73 

8.65±5.48 

4.14+4.70 

5.66+2.94 

0.69  0.03 

0.17 

Uterine  weight  (g) 

0.10±0.05 

0.04±0.02 

0.10+0.02 

0.05+0.02 

0.30  <0.0001 

0.34 

Distal  Femur 

BV/TV 

0.019±0.004 

0.012±0.004 

0.038+0.011 

0.021+0.010 

<0.0001  <0.0001 

0.06 

SMI 

3.818±0.250 

4.055±0.357 

3.387+0.263 

3.810+0.294 

0.0008  0.0011 

0.32 

Tb.N  (mm-') 

2.554±0.239 

2.165±0.385 

3.128+0.218 

2.797+0.276 

<0.0001  0.0004 

0.76 

Tb.Th  (mm) 

0.040±0.005 

0.041±0.005 

0.039+0.003 

0.037+0.004 

0.08  0.87 

0.23 

Tb.Sp  (mm) 

0.393±0.036 

0.477±0.097 

0.317+0.026 

0.359+0.037 

<0.0001  0.0012 

0.25 

L5  Vertebra 

BV/TV 

0.189±0.028 

0.177±0.013 

0.253+0.019 

0.212+0.019 

<0.0001  0.0004 

0.05 

Tb.N  (mm-') 

3.797±0.513 

3.580±0.285 

4.491+0.345 

4.022+0.254 

<0.0001  0.0091 

0.32 

Tb.Th  (mm) 

0.050±0.003 

049±0.002 

0.056+0.002 

0.053+0.002 

<0.0001  0.0032 

0.02 

Tb.Sp  (mm) 

0.227±0.023 

0.229±0.013 

0.202+0.020 

0.214+0.012 

0.0013  0.25 

0.44 

Serum  PINP 

WT  :  Amp4  ^ 

n-values  IG;  T;  G  x  TU^^^ 

n-values  IG;  T;  G  x  Tl'^'"'^^ 

Pre-op 

Post- 

,Op4wks 

Post-op'^"'^® 

5.62±1.21  :  6.02±1.41 

5.83±1.41  :  4.77±1.22  3.15+0.65  :  2.81±0.76 

0.34;  0.16;  0.05 

0.95;  <0.0001;  0.36 

Serum  CTX 

WT  :  Nmp4  '^ 

n-values  IG;  T;  G  x  TU^^^ 

n-values  IG;  T;  G  x  Tl'^'"’^^ 

Pre-op 

Post- 

Op4wks 

Post-op'^"'"® 

13.66±2.43  :  13.37±1.88 

12.96+3.04  :  13.53+2.58  11.47+2.24  :  9.36+1.22 

0.92;  0.46;  0.96 

0.14;  0.0003;  0.18 

The  %  change  body  weight,  uterine  weight,  and  microCT  of  distal  femur  and  L5  vertebra  from  WT  and  Nmp4'^^  mice  after  ovx  or  sham  operation  4wks  post-op. 
Serum  bone  formation  [PINP]  and  bone  resorption  [CTX]  marker  levels  were  compared  in  mice  previous  to  the  operation  (Pre-op)  and  4wks  and  12wks 
subsequent  to  surgery  (Post-op).  The  12wks  post-op  data  was  obtained  from  the  vehicle-control  treatment  groups.  Data  are  average  ±  SD,  number  of 
mice/experimental  group=8-14  [4  mice  in  WT  SHAM  uterine  weight]).  Statistical  significance  was  set  at  p<0.05  and  differences  were  determined  using  a  2-way 
ANOVA. 
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Table  2:  PTH-induced  bone  gain  data. 


WT 

Nmp4^^ 

2-WAY  ANOVA  n-values 

Distal  Femur 

4wks 

VEH 

PTH 

VEH 

PTH 

Genotype 

Treatment 

Gene  x  Treat 

Conn  D  (mm'^) 

3.180±3.870 

33.230±26.730 

9.681±15.979 

67.533±14.111 

0.0018 

<0.0001 

0.03 

SMI 

3.752±0.437 

3.013±0.384 

3.472±0.327 

2.514±0.113 

0.0025 

<0.0001 

0.36 

Tb.N  (mm-') 

2.100±0.519 

2.441±0.281 

2.712±0.241 

2.833±0.224 

0.0002 

0.06 

0.36 

Tb.Th  (mm) 

0.039±0.010 

0.042±0.007 

0.033±0.003 

0.044±0.003 

0.54 

0.004 

0.09 

Tb.Sp  (mm) 

8wks 

0.510±0.157 

0.409±0.051 

0.370±0.036 

0.342±0.032 

0.0019 

0.04 

0.24 

Conn  D  (mm"^) 

3.123±5.307 

38.658±14.910 

0.982±1.103 

58.128±13.570 

0.03 

<0.0001 

0.0064 

SMI 

3.808±0.479 

2.470±0.284 

3.589±0.218 

2.262±0.141 

0.05 

<0.0001 

0.96 

Tb.N  (mm-') 

2.132±0.297 

2.164±0.431 

2.286±0.145 

2.552±0.277 

0.02 

0.17 

0.28 

Tb.Th  (mm) 

0.037±0.006 

0.048±0.005 

0.030±0.004 

0.049±0.003 

0.12 

<0.0001 

0.02 

Tb.Sp  (mm) 

0.476±0.072 

0.471±0.109 

0.438±0.033 

0.378±0.045 

0.01 

0.20 

0.27 

L5  Vertebra 

4wks 

VEH 

PTH 

VEH 

PTH 

Genotype 

Treatment 

Gene  x  Treat 

Tb.N  (mm-') 

3.453±0.451 

4.875±0.587 

3.891±0.504 

5.518±0.381 

0.0049 

<0.0001 

0.56 

Tb.Th  (mm) 

0.051±0.002 

0.049±0.002 

0.054±0.004 

O.OSliO.OOl 

0.04 

0.03 

0.60 

Tb.Sp  (mm) 

8wks 

0.246±0.021 

0.224±0.030 

0.229±0.021 

0.197±0.021 

0.02 

0.0036 

0.52 

Tb.N  (mm-') 

4.046±0.917 

5.648±1.191 

3.627±0.235 

5.906±0.754 

0.79 

<0.0001 

0.26 

Tb.Th  (mm) 

0.053±0.003 

0.049±0.004 

0.055±0.001 

0.054±0.001 

0.0018 

0.0044 

0.09 

Tb.Sp  (mm) 

0.239±0.021 

0.206±0.037 

0.256±0.020 

0.186±0.023 

0.86 

<0.0001 

0.05 

MicroCT  (distal  femur  and  L5  vertebra)  from  ovx  WT  and  ovx  Nmp4  '^  mice  after  4  wks  and  8wks  PTH/VEH  therapy.  Data  are  average  ±  SD,  number  of 
mice/experimental  group  =  8-9.  Statistical  significance  was  set  at  p<0.05  and  differences  were  determined  using  a  2-way  ANOVA. 
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Table  3:  Histomorphometry  and  serum  analyses. 


WT _  Nmp4-'- _  _ 2-WAY  ANOVA  p-values 


Dynamic  histo 

VEH 

PTH 

VEH 

PTH 

Genotype 

Treatment 

Gene  x  Treat 

MAR  (pm/day) 

2.28±0.37 

3.80±0.73 

2.29±0.37 

3.61±0.40 

0.70 

<0.0001 

0.66 

MS/BS  (%) 

0.41±0.09 

0.55±0.05 

0.44±0.10 

0.52±0.06 

0.98 

0.01 

0.45 

BFR  (pmVpm/day) 

0.95+0.28 

2.09+0.52 

1.01+0.25 

1.86+0.22 

0.60 

<0.0001 

0.37 

Serum 


PINP  (ng/ml) 
CTX  (ng/ml) 


WT 

Nmp4'^^ 

VEH^^ks 

pTpjSwks 

VEH^^ks 

pTpjSwks 

Genotype 

3.147+0.653 

10.066F2.659 

2.806+0.760 

8.042F3.304 

0.19 

11.466±2.239 

15.147+3.518 

9.361F1.222 

14.157±1.532 

0.12 

2-WAY  ANOVA  p-values _ 

Treatment  Gene  x  Treat 

<0.0001  0.34 

0.0002  0.56 


Dynamic  bone  histomorphometry  data  of  the  distal  femur  from  WT  and  Nmp4''^  mice  treated  with  intermittent  PTH  or  vehicle  for  4wks  (8wks  post-op).  Sera  data 
were  collected  at  the  end  of  8wks  treatment  (12wks  post-op).  The  parameters  include  mineral  apposition  rate  (MAR),  mineralizing  surface/hone  surface  (MS/BS), 
and  hone  formation  rate  (BFR).  Data  are  average  +  SD,  number  of  mice/experimental  group  =  4-7.  A  2-way  ANOVA  was  used  to  determine  statistical  differences 
and  significance  was  set  at  p<0.05. 
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Table  4:  ENCODE  ChIP-Seq  Significance  Tool  profile  for  enriched  transcription  factors  [TFs]  within  the  Nmp4  target  core  gene  list 


Factor 

0-vaIue* 

Factor 

0-vaIue 

Nmp4 

O.OOE+00 

Max 

O.OOE+00 

CHD2 

O.OOE+00 

Mxil 

O.OOE+00 

CTCF 

O.OOE+00 

NELFe 

O.OOE+00 

GCN5 

O.OOE+00 

Pol2 

O.OOE+00 

HCFCl 

O.OOE+00 

SIN3A 

O.OOE+00 

MAZ 

O.OOE+00 

TBP 

O.OOE+00 

p300 

O.OOE+00 

c-Myc 

7.352e-317 

*  Hypergeometric  test;  Benjamini-Hochberg;  (select  TFs  from  72  entries). 
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Table  5:  DAVID  profile  of  KEGG  pathway  mapping. 


GO  Term  Pathways _ FDR 

TOR  signaling  pathway  0.003 

Insulin  signaling  pathway  0.004 

Chronic  myeloid  leukemia  0.026 

JAK-STAT  signaling  pathway  0.026 

Neurotrophin  signaling  pathway  0.034 


Only  pathways  with  an  FDR  of  p<0.05  are  listed 
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1  FIGURE  LEGENDS 

2 

3  Figure  1:  Schematic  of  treatment  regimen  for  WT  and  Nmp4^^^  mice;  Group  1  mice  were  subjected  to 

4  ovariectomy  (ovx)  or  sham  operation  at  12wks  of  age  and  evaluated  for  bone  loss  4wks  post-op  (16wks  of 

5  age).  Group  2  mice  were  ovx  at  12wks  of  age  and  began  PTH  or  vehicle  therapy  at  16wks  of  age  for  a 

6  duration  of  4wks  and  8wks.  Endpoint  analyses  included  micro-computed  tomography  pCT,  serum 

7  analysis  for  N-terminal  propeptide  of  type  1  procollagen  (PINP)  and  C-terminal  telopeptides  (CTX),  and 

8  dynamic  histomorphometry. 

9 

10  Figure  2:  Disabling  Nmp4  enhances  PTH  restorative  therapy  in  the  distal  femur  of  ovx  Nmp4^^^  mice  [A] 

11  Interaction  plots  of  femoral  trabecular  bone  volume/total  volume  (BV/TV)  of  ovx  WT  and  ovx  Nmp4^^^ 

12  mice  as  determined  by  pCT  at  4wks  of  treatment  and  8wks  of  treatment.  Data  are  average  +  SD,  number 

13  of  mice/experimental  group  =  8-9).  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

14  significance  was  set  at  p<0.05.  The  Tukey’s  HSD  post  hoc  test  was  used  to  determine  differences  between 

15  the  treatment  groups.  There  were  genotype,  treatment  and  genotype  x  treatment  interaction  at  both  time 

16  points.  There  was  no  difference  between  the  vehicle-treated  WT  and  Amp4  ^  mice.  [B]  pCT  images 

1 7  showing  PTH-induced  improvements  in  distal  femur  trabecular  architecture  in  ovx  WT  and  Amp4  ^  mice 

18  after  8  weeks  of  treatment  (12wks  post-op,  24wks  of  age). 

19 

20  Figure  3:  The  exaggerated  response  to  anabolic  PTH  persists  in  the  L5  vertebra  of  ovx  Nmp4^'^  mice.  [A] 

21  Interaction  plots  of  L5  vertebra  bone  volume/total  volume  (BV/TV)  of  ovx  WT  and  ovx  ^  mice  as 

22  determined  by  pCT  at  4wks  of  treatment  and  8wks  of  treatment.  Data  are  average  +  SD,  number  of 

23  mice/experimental  group  =  8-9).  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

24  significance  was  set  at  p<0.05.  The  LS  Means  Student  t  post  hoc  test  was  used  to  determine  differences 

25  between  the  treatment  groups.  There  were  genotype,  treatment  effects  at  both  time  points  and  a  genotype 

26  X  treatment  interaction  at  8wks  therapy.  There  was  no  difference  between  the  vehicle-treated  WT  and 
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27  Nmp4'^^  mice.  [B]  |j,CT  images  showing  PTH-induced  improvements  in  L5  trabecular  architecture  in  ovx 

28  WT  and  ^  mice  after  8  weeks  of  treatment  (12wks  post-op,  24wks  of  age). 

29 

30  Figure  4  :  Ovx  does  not  abrogate  the  expanded  population  of  osteoprogenitors  and  CD8-I-  T  cells  in  Nmp4' 

31  ^  mice.  FACS  analysis  of  BM  and  PBL  osteoprogenitors,  CD8-I-  T  cells,  and  CD4-I-  T  cells.  [A,  D]  The 

32  frequency  of  femoral  BM  and  PBL  CD45-/CD105-i-/CD146-i-/CD105-i-/nestin-i-  osteoprogenitor  cells  in 

33  WT  and  Amp4  ^  mice  at  the  end  of  4wks  and  8wks  treatment  with  intermittent  PTH  or  vehicle  control;  [B, 

34  E]  the  frequency  of  BM  and  PBL  CD8-I-  T  cells  from  the  WT  and  Amp4  ^  mice;  [C,  F]  the  frequency  of 

35  BM  and  PBL  CD4-I-  T  cells  from  the  WT  and  null  mice.  Data  are  average  +  SD,  number  of 

36  mice/experimental  group  =  8-9;  Statistical  differences  were  determined  using  a  2-way  ANOVA  and 

37  significance  was  set  at  p<0.05. 

38 

39  Figure  5:  Expanded  Amp4  ^  MSPCs  exhibit  enhanced  proliferation  and  mineralization  in  culture.  [A] 

40  Comparative  growth  rates  of  expanded  WT  and  Amp4  ^  MSPCs.  Cell  counts/day  (n=4  lines  per  genotype 

41  loglO  cells/well,  3  wells/sample,  average  +  SD,  t  test,  t<0.05).  Note:  each  ‘line’  is  derived  from  a  single 

42  mouse  [B]  Alkaline  phosphatase  (alk  phos)  and  alizarin  red  staining  of  a  WT  and  Amp4  ^  MSPC  cultures 

43  from  Day7-Day28.  See  text  for  details 

44 

45  Figure  6:  Nmp4  associates  with  core  target  genes  common  to  multiple  cell  types  and  acts  as  a  negative 

46  regulation  of  cellular  biosynthetic  processes  [A]  Venn  diagram  illustrating  the  shared  Nmp4  target  genes 

47  in  the  MC3T3-E1  osteoblast-like  cells  (vehicle-treated),  and  the  three  ENCODE  cells  lines,  ES-E14 

48  (embryonic  stem  cells),  MEL,  and  CH12  cells  (B-cell  lymphomas).  [B]  DAVID/REVIGO  gene  ontology 

49  (GO)  profile  of  Nmp4  core  target  genes 

50 

51  Figure  7:  Nmp4  binds  to  AT-rich  DNA  typically  proximal  to  TSS  sites  or  within  intragenic  regions.  [A] 

52  Genome-wide  mapping  of  the  Nmp4  binding  sites  show  that  most  sites  are  distributed  in  the  TSS  and 
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intragenic  regions  of  the  genome.  ChIP-seq  analysis  included  vehicle-treated  and  PTH-treated  MC3T3-E1 
osteoblast-like  cells  (vMC  and  pMC,  respectively)  and  three  murine  cell  lines  from  the  ENCODE 
Consortium  including  ES-E14  (Esl4),  which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and 
two  mouse  erythroleukemia  cell  lines  (Chl2  and  MEL)  derived  from  B-cell  lymphomas.  [B]  GEM 
analysis  for  the  Nmp4  consensus  sequence  derived  from  MC3T3-E1  cells.  A  minimal  k-mer  width  of  6 
and  maximum  of  20  were  used.  The  optimal  position  weight  matrix  (PWM)  score  for  the  MC3T3-E1  data 
was  10.07.  The  hypergeometric  P-value  (hgp)  was  le-1466.1. 

Figure  8:  ChIP-seq  reveals  Nmp4  binding  profiles  at  specific  gene  loci.  Mouse  MC3T3-E1  cells  were 
seeded  into  twenty-one  150mm  plates  at  an  initial  density  of  50,000  cells/plate  (320  cells/cm^)  and 
maintained  in  aMEM  complete  medium  -i-  ascorbic  acid  for  14  days.  Prior  to  harvest  cells  were  treated 
with  25nM  hPTH(l-34)  or  vehicle  control  for  Ihr.  Processing  for  ChIP-seq  analysis  was  performed  as 
described  in  the  Materials  and  Methods.  Sequences  (50nt  reads,  single  end)  were  aligned  to  the  mouse 
genome  (mmlO)  using  the  BWA  algorithm.  Alignments  were  extended  in  silico  at  their  3’-ends  to  a 
length  of  150bp,  which  is  the  average  genomic  fragment  length  in  the  size-selected  library,  and  assigned 
to  32-nt  bins  along  the  genome.  Nmp4  (Znf384)  peak  locations  were  determined  using  the  MACS 
algorithm  (vl.4.2)  with  a  cutoff  of  pvalue  =  le-7.  The  genomic  loci  including  the  chromosome  number 
and  nucleotide  interval  are  indicated.  Read  scales  are  indicated  on  the  Y-axis.  An  arrow  indicates  the 
transcriptional  start  sites  and  direction  of  transcription  for  each  of  the  genes;  vertical  boxes  within  the 
gene  indicate  exons.  The  Nmp4  ChIP-seq  gene  profiles  include  (A)  Nid2  (B)  Akt2,  (C)  Pdkl  (D)  ccdc53, 
(E)  Arrb2  and  (E)  Irsl.  The  input  DNA  profiles  were  devoid  of  peaks. 

Figure  9:  Comparison  of  mRNA  expression  profiles  derived  from  non-differentiating  (Day  3)  and 
osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^'  cells.  All  transcript  levels  are  compared  to  WT 
Day  3  providing  a  time  course  of  expression.  mRNA  profiles  [A]  Igfbp2-,  [B]  lgfbp4;  [C]  Pdkl  ',  [D] 

Bmp2\  [E]  Pthlr  were  derived  from  the  TED  A  system  (Eormat  96a,  Applied  Biosystems,  Eoster  City, 
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CA)  performed  on  a  QuantStudio™  7  Flex  Real-Time  PCR  System  and  normalized  with  GusB.  Profile 
[F]  Bglap  mRNA  profile  qRT-PCR  reactions  were  performed  on  an  Eppendorf  Mastercycler®  RealPlex^ 
using  Rplpl  Mm03059047_gH)  as  the  normalizer  as  previously  described  [Robling  et  al,  2009]. 
Comparison  of  profiles  using  GusB  and  Rplpl  as  the  normalizer  showed  no  differences  in  the  shape  of  the 
expression  profiles 

Figure  10:  Comparison  of  mRNA  expression  profiles  derived  from  non-differentiating  (Day  3)  and 
osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^~  cells.  All  transcript  levels  are  compared  to  WT 
Day  3  providing  a  time  course  of  expression.  mRNA  profiles  [A]  Cxcll2\  [B]  Plaur;  [C]  Sppl',  [D] 

Thbsl',  [E]  Collal  were  derived  from  the  TLDA  system  (Eormat  96a,  Applied  Biosystems,  Eoster  City, 
CA)  performed  on  a  QuantStudio™  7  Elex  Real-Time  PCR  System  and  normalized  with  GusB.  Profile 
[E]  Sp7  mRNA  profile  qRT-PCR  reactions  were  performed  on  an  Eppendorf  Mastercycler®  RealPlex^ 
using  Rplpl  (Mm03059047_gH)  as  the  normalizer  as  previously  described  [Robling  et  al.,  2009].  The 
Day  16  WT  sample  is  the  average  of  two  replicates.  Comparison  of  profiles  using  GusB  and  Rplp2  as  the 
normalizer  showed  no  differences  in  the  shape  of  the  expression  profiles 

Supplemental  Figure  SI:  qRT-PCR  validates  the  ChIP-seq  profiles.  [A]  The  Nmp4  ChIP-seq  profile  for 
the  gene  Collal.  The  genomic  loci  including  the  chromosome  number  and  nucleotide  interval  are 
indicated.  Read  scale  is  indicated  on  the  Y-axis.  An  arrow  marks  the  transcriptional  start  site  and  direction 
of  transcription;  vertical  boxes  within  the  gene  identify  exons.  [B]  qRT-PCR  was  used  to  authenticate  the 
ChIP-seq  peaks  as  described  in  the  Materials  and  Methods. 

Supplemental  Figure  2:  Volcano  plots  derived  from  gene  expression  profiles  of  non-differentiating  (Day 
3)  and  osteogenic-differentiating  (Days  5-16)  WT  and  Nmp4'^'  as  described  in  the  Materials  and  Methods. 
RNA  expression  profiling  was  performed  on  a  QuantStudio™  7  Flex  Real-Time  PCR  System  and  data 
analyzed  using  the  ExpressionSuite  v  1.0.4™  analysis  software  (Applied  Biosystems)  as  described  in  the 
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Materials  and  Methods.  [A]  WT  vs  Nmp4'^'  cells  at  Day  3  post-seeding.  mRNA  transcript  expression  was 
compared  to  WT  cells  (Day  3).  Cells  maintained  in  Mesencult™  Media  -i-  Mesencult™  Stimulatory 
Supplement.  [B]  WT  vs  Nmp4'^'  cells  at  Day  5  post-seeding.  mRNA  transcript  expression  was  compared 
to  WT  cells  (Day  5).  Cells  maintained  in  differentiation  medium  for  48hrs  [C]  WT  vs  Nmp4'^^  cells  at  Day 
7  post-seeding.  mRNA  transcript  expression  was  compared  to  WT  cells  (Day  7).  Cells  maintained  in 
differentiation  medium  for  96hrs  [D]  WT  vs  Nmp4'^'  cells  at  Day  9  post-seeding.  mRNA  transcript 
expression  was  compared  to  WT  cells  (Day  9).  Cells  maintained  in  differentiation  medium  for  144hrs  [E] 
WT  vs  Nmp4~^'  cells  at  Day  16  post-seeding.  mRNA  transcript  expression  was  compared  to  WT  cells 
(Day  16).  Cells  maintained  in  differentiation  medium  for  192hrs.  Genes  indicated  with  green  dots  (left  of 
center)  and  above  the  X-axis  exhibited  a  significant  downregulation  by  over  2-fold.  Genes  indicated  with 
the  red  dots  (right  of  center)  and  above  the  X-axis  exhibited  a  significant  upregulation  by  over  2-fold. 

Supplemental  TABLE  1  :  96  Nmp4  ‘core’  target  genes,  non-core  target  genes,  and  non-target  genes 
including  5  candidate  normalizer  genes.  Individual  cDNAs  were  quantified  by  qRT-PCR  using  a  custom 
TLDA  system  (Format  96a,  Applied  Biosystems,  Foster  City,  CA)  as  described  in  the  Materials  and 
Methods  on  a  QuantStudio™  7  Flex  Real-Time  PCR  System.  We  used  the  ExpressionSuite  vl.0.4™ 
analysis  software  (Applied  Biosystems)  to  analyze  these  data. 

Supplemental  TABLE  2  :  21 14  Nmp4  ‘core’  target  genes  common  to  the  four  cell  lines  MC3T3-E1 
osteoblast-like  cells,  and  three  murine  cell  lines  from  the  ENCODE  Consortium  including  ES-E14  (Esl4), 
which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and  two  mouse  erythroleukemia  cell  lines 
(Chi 2  and  MEL)  derived  from  B-cell  lymphomas. 

Supplemental  TABLE  3  :  IPA  analysis  of  21 14  Nmp4  ‘core’  target  genes  common  to  the  four  cell  lines 
MC3T3-E1  osteoblast-like  cells,  and  three  murine  cell  lines  from  the  ENCODE  Consortium  including 


5 


130  ES-E14  (Esl4),  which  are  E14  undifferentiated  mouse  embryoinic  stem  cells,  and  two  mouse 

131  erythroleukemia  cell  lines  (Chl2  and  MEL)  derived  from  B-cell  lymphomas. 
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Detector 

Reporter 

Quencher  Description 

Comments  Sequence 

18S-Hs99999901_sl 

FAM 

Non  Fluorescent 

Eukaryotic  18S  rRNA 

Akt2-Mm02026778_gl 

FAM 

Non  Fluorescent 

thymoma  viral  proto-oncogene  2 

Alpl-Mm00475834_ml 

FAM 

Non  Fluorescent 

alkaline  phosphatase,  liver/bone/kidney 

Arrb2-Mm00520666_gl 

FAM 

Non  Fluorescent 

arrestin,  beta  2 

Atf4-Mm00515325_gl 

FAM 

Non  Fluorescent 

activating  transcription  factor  4 

Bcl2lll-Mm00437796_ml 

FAM 

Non  Fluorescent 

BCL2-like  11  (apoptosis  facilitator) 

Bglap3-Mm00649782_gH 

FAM 

Non  Fluorescent 

bone  gamma-carboxyglutamate  protein  3 

Bmp2-Mm01340178_ml 

FAM 

Non  Fluorescent 

bone  morphogenetic  protein  2 

Bmp3-Mm00557790_ml 

FAM 

Non  Fluorescent 

bone  morphogenetic  protein  3 

Bmp4-Mm00432087_ml 

FAM 

Non  Fluorescent 

bone  morphogenetic  protein  4 

Bmp6-Mm01332882_ml 

FAM 

Non  Fluorescent 

bone  morphogenetic  protein  6 

Bmprla-Mm00477650_ml 

FAM 

Non  Fluorescent 

bone  morphogenetic  protein  receptor,  type  lA 

Cbfb-Mm01251026_gl 

FAM 

Non  Fluorescent 

core  binding  factor  beta 

Cbx4-Mm00483089_ml 

FAM 

Non  Fluorescent 

chromobox  4 

Cbx7-Mm00520006_m  1 

FAM 

Non  Fluorescent 

chromobox  7 

Cdc6-Mm03048221_ml 

FAM 

Non  Fluorescent 

cell  division  cycle  6 

Cirbp-Mm00483336_gl 

FAM 

Non  Fluorescent 

cold  inducible  RNA  binding  protein 

Collal-Mm00801666_gl 

FAM 

Non  Fluorescent 

collagen,  type  1,  alpha  1 

Crebbp-Mm01342452_ml 

FAM 

Non  Fluorescent 

CREB  binding  protein 

Cxcll2-Mm00445553_ml 

FAM 

Non  Fluorescent 

chemokine  (C-X-C  motif)  ligand  12 

Ddit4-Mm00512504_gl 

FAM 

Non  Fluorescent 

DNA-damage-inducible  transcript  4 

Dnm3os-Mm03455916_sl 

FAM 

Non  Fluorescent 

dynamin  3,  opposite  strand 

Ehmt2-Mm01132261_ml 

FAM 

Non  Fluorescent 

euchromatic  histone  lysine  N-methyltransferase  2 

Eif4e-Mm00725633_sl 

FAM 

Non  Fluorescent 

eukaryotic  translation  initiation  factor  4E 

Eif5a-Mm01971736_gl 

FAM 

Non  Fluorescent 

eukaryotic  translation  initiation  factor  5A 

Ephb4-Mm01201157_ml 

FAM 

Non  Fluorescent 

Eph  receptor  B4 

Ezhl-Mm00468440_ml 

FAM 

Non  Fluorescent 

enhancer  of  zeste  homolog  1  (Drosophila) 

Fos-Mm00487425_ml 

FAM 

Non  Fluorescent 

FBJ  osteosarcoma  oncogene 

Fosl2-Mm00484442_ml 

FAM 

Non  Fluorescent 

fos-like  antigen  2 

Gadd45b-Mm00435121_gl 

FAM 

Non  Fluorescent 

growth  arrest  and  DNA-damage-inducible  45  beta 

Gadd45g-Mm01352550_gl 

FAM 

Non  Fluorescent 

growth  arrest  and  DNA-damage-inducible  45  gamma 

Gasl-Mm01700206_gl 

FAM 

Non  Fluorescent 

growth  arrest  specific  1 

Hdac3-Mm00515916_ml 

FAM 

Non  Fluorescent 

histone  deacetylase  3 

Hifla-Mm00468869_ml 

FAM 

Non  Fluorescent 

Igfl-Mm00439560_ml 

FAM 

Non  Fluorescent 

Igflr-Mm00802831_ml 

FAM 

Non  Fluorescent 

Igfbp2-Mm00492632_ml 

FAM 

Non  Fluorescent 

Igfbp4-Mm00494922_ml 

FAM 

Non  Fluorescent 

Igfbp6-Mm00599696_ml 

FAM 

Non  Fluorescent 

Illrl-Mm00434237_ml 

FAM 

Non  Fluorescent 

Ing3-Mm00458324_ml 

FAM 

Non  Fluorescent 

Ing4-Mm00460097_ml 

FAM 

Non  Fluorescent 

Insr-Mm01211875_ml 

FAM 

Non  Fluorescent 

Irsl-Mm01278327_ml 

FAM 

Non  Fluorescent 

Klf4-Mm00516104_ml 

FAM 

Non  Fluorescent 

Kpna2-Mm00834020_gH 

FAM 

Non  Fluorescent 

Lgals3-Mm00802901_ml 

FAM 

Non  Fluorescent 

Lrpl-Mm00464608_ml 

FAM 

Non  Fluorescent 

Msx2-Mm00442992_ml 

FAM 

Non  Fluorescent 

Myc-Mm00487804_ml 

FAM 

Non  Fluorescent 

Nbrl-Mm01249798_ml 

FAM 

Non  Fluorescent 

Ncorl-Mm01333102_ml 

FAM 

Non  Fluorescent 

Ncor2-Mm00448796_ml 

FAM 

Non  Fluorescent 

Nfil3-Mm00600292_sl 

FAM 

Non  Fluorescent 

Nr3cl-Mm00433832_ml 

FAM 

Non  Fluorescent 

Pdkl-Mm00554306_ml 

FAM 

Non  Fluorescent 

Perl-Mm00501813_ml 

FAM 

Non  Fluorescent 

Phfl2-Mm00663497_ml 

FAM 

Non  Fluorescent 

Piml-Mm00435712_ml 

FAM 

Non  Fluorescent 

Plaur-Mm00440911_ml 

FAM 

Non  Fluorescent 

Pou2fl-Mm00448332_ml 

FAM 

Non  Fluorescent 

Ppplrl5a-Mm01205601_gl 

FAM 

Non  Fluorescent 

Prpf  19-M  m00467298_m  1 

FAM 

Non  Fluorescent 

Psmal-Mm00803741_ml 

FAM 

Non  Fluorescent 

Psma2-Mm00776364_mH 

FAM 

Non  Fluorescent 

Psma3-Mm00834115_gl 

FAM 

Non  Fluorescent 

Psmb6-Mm00599713_gl 

FAM 

Non  Fluorescent 

hypoxia  inducible  factor  1,  alpha  subunit 

insulin-like  growth  factor  1 

insulin-like  growth  factor  I  receptor 

insulin-like  growth  factor  binding  protein  2 

insulin-like  growth  factor  binding  protein  4 

insulin-like  growth  factor  binding  protein  6 

interleukin  1  receptor,  type  I 

inhibitor  of  growth  family,  member  3 

inhibitor  of  growth  family,  member  4 

insulin  receptor 

insulin  receptor  substrate  1 

Kruppel-like  factor  4  (gut) 

karyopherin  (importin)  alpha  2 

lectin,  galactose  binding,  soluble  3 

low  density  lipoprotein  receptor-related  protein  1 

msh  homeobox  2 

myelocytomatosis  oncogene 

neighbor  of  Brcal  gene  1 

nuclear  receptor  co-repressor  1 

nuclear  receptor  co-repressor  2 

nuclear  factor,  interleukin  3,  regulated 

nuclear  receptor  subfamily  3,  group  C,  member  1 

pyruvate  dehydrogenase  kinase,  isoenzyme  1 

period  circadian  clock  1 

PHD  finger  protein  12 

proviral  integration  site  1 

plasminogen  activator,  urokinase  receptor 

POU  domain,  class  2,  transcription  factor  1 

protein  phosphatase  1,  regulatory  (inhibitor)  subunit  15A 

PRP19/PS04  pre-mRNA  processing  factor  19  homolog  (S.  cerevisiae) 

proteasome  (prosome,  macropain)  subunit,  alpha  type  1 

proteasome  (prosome,  macropain)  subunit,  alpha  type  2 

proteasome  (prosome,  macropain)  subunit,  alpha  type  3 

proteasome  (prosome,  macropain)  subunit,  beta  type  6 


Pthlr-Mm00441046_ml 

FAM 

Non  Fluorescent 

Ptp4al-Mm00850755_gl 

FAM 

Non  Fluorescent 

Rad23b-Mm00772280_ml 

FAM 

Non  Fluorescent 

Rffl-Mm00482724_ml 

FAM 

Non  Fluorescent 

Rnfl67-Mm00550967_gl 

FAM 

Non  Fluorescent 

Rnf5-Mm01134793_gl 

FAM 

Non  Fluorescent 

Rps6kbl-Mm01310033_ml 

FAM 

Non  Fluorescent 

Rptor-Mm00712676_ml 

FAM 

Non  Fluorescent 

Runx2-Mm00501584_ml 

FAM 

Non  Fluorescent 

Sin3a-Mm00488255_ml 

FAM 

Non  Fluorescent 

Smad7-Mm00484742_ml 

FAM 

Non  Fluorescent 

Sp7-Mm04209856_ml 

FAM 

Non  Fluorescent 

Sparc-Mm00486332_ml 

FAM 

Non  Fluorescent 

Sppl-Mm00436767_ml 

FAM 

Non  Fluorescent 

Suzl2-Mm01304145_gl 

FAM 

Non  Fluorescent 

Thbs2-Mm01279240_m  1 

FAM 

Non  Fluorescent 

Tle3-Mm00437097_ml 

FAM 

Non  Fluorescent 

Tob2-Mm00451524_sl 

FAM 

Non  Fluorescent 

Ttc3-Mm00493917_ml 

FAM 

Non  Fluorescent 

Ubr2-Mm00524868_ml 

FAM 

Non  Fluorescent 

Uspl5-Mm00452856_ml 

FAM 

Non  Fluorescent 

Usp2-Mm00497452_ml 

FAM 

Non  Fluorescent 

Xiap-Mm01311594_mH 

FAM 

Non  Fluorescent 

Zbtb7a-Mm00657132_ml 

FAM 

Non  Fluorescent 

Rplp2-Mm00782638_sl 

FAM 

Non  Fluorescent 

Gusb-Mm01197698_ml 

FAM 

Non  Fluorescent 

B2m-Mm00437762_ml 

FAM 

Non  Fluorescent 

Hprt-Mm01545399_ml 

FAM 

Non  Fluorescent 

Dkk2-Mm01322146  ml 

FAM 

Non  Fluorescent 

parathyroid  hormone  1  receptor 

protein  tyrosine  phosphatase  4al 

RAD23b  homolog  (S.  cerevisiae) 

ring  finger  and  FYVE  like  domain  containing  protein 

ring  finger  protein  167 

ring  finger  protein  5 

ribosomal  protein  S6  kinase,  polypeptide  1 
regulatory  associated  protein  of  MTOR,  complex  1 
runt  related  transcription  factor  2 
transcriptional  regulator,  SIN3A  (yeast) 

SMAD  family  member  7 

Sp7  transcription  factor  7 

secreted  acidic  cysteine  rich  glycoprotein 

secreted  phosphoprotein  1 

suppressor  of  zeste  12  homolog  (Drosophila) 

thrombospondin  2 

transducin-like  enhancer  of  split  3,  homolog  of  Drosophila  E(spl) 

transducer  of  ERBB2,  2 

tetratricopeptide  repeat  domain  3 

ubiquitin  protein  ligase  E3  component  n-recognin  2 

ubiquitin  specific  peptidase  15 

ubiquitin  specific  peptidase  2 

X-linked  inhibitor  of  apoptosis 

zinc  finger  and  BTB  domain  containing  7a 

ribosomal  protein,  large  P2 

glucuronidase,  beta 

beta-2  microglobulin 

hypoxanthine  guanine  phosphoribosyl  transferase 
dickkopf  homolog  2  (Xenopus  laevis) 
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Common  elements  in  "vMC3T3",  "MEL",  "Esl4"  and  "Chl2":  Nmp4  'CORE'  genes 

0610031016Rik 

0610043K17Rik 

1110001J03Rik 

1110002J07Rik 

1110002L01Rik 

1110037F02Rik 

1110038B12Rik 

1110051M20Rik 

1190002F15Rik 

1600002K03Rik 

1600020E01Rik 

1700001GllRik 

1700007K13Rik 

1700007L15Rik 

1700013F07Rik 

1700016C15Rik 

1700018L02Rik 

1700018M17Rik 

1700021A07Rik 

1700022K14Rik 

1700023H06Rik 

1700030C12Rik 

1700034H15Rik 

1700052N19Rik 

1700063D05Rik 

1700064E03Rik 

1700067K01Rik 

1700095J07Rik 


1700101lllRik 

1700110C19Rik 

1700112E06Rik 

1700120B22Rik 

1810013L24Rik 

1810019N24Rik 

1810059C17Rik 

2010204K13Rik 

2210016F16Rik 

2210408F21Rik 

2210417K05Rik 

2310011J03Rik 

2310034G01Rik 

2310068J16Rik 

2410002F23Rik 

2410022MllRik 

2500002B13Rik 

2500004C02Rik 

2510009E07Rik 

2610301B20Rik 

2610507BllRik 

2700029M09Rik 

2700050L05Rik 

2700060E02Rik 

2810403A07Rik 

2810403D21Rik 

2810404F17Rik 

2810408M09Rik 

2810428ll5Rik 


2810454H06Rik 

3010003L21Rik 

3110009E18Rik 

3110067C02Rik 

3110082J24Rik 

4732471J01Rik 

4921530L18Rik 

4921531C22Rik 

4930402F06Rik 

4930404l05Rik 

4930412O13Rik 

4930432B10Rik 

4930447C04Rik 

4930503L19Rik 

4930509E16Rik 

4930529M08Rik 

4930546H06Rik 

4930552N02Rik 

4931406C07Rik 

4933411K20Rik 

4933433G15Rik 

4933434E20Rik 

5031434011Rik 

5330430P22Rik 

5530601H04Rik 

5730405O15Rik 

5730420D15Rik 

5730455P16Rik 

5730508B09Rik 


9030624J02Rik 

9430008C03Rik 

9430041J12Rik 

9430083A17Rik 

9530026F06Rik 

9530027J09Rik 

9530068E07Rik 

9630014M24Rik 

A330017A19Rik 

A330050B17Rik 

A330069E16Rik 

A430018G15Rik 

A630072M18Rik 

A730036ll7Rik 

A830031A19Rik 

A830035A12Rik 

A930001C03Rik 

A930007ll9Rik 

AA387883 

AA415398 

AA465934 

Aaas 

Aarsdl 

AB041803 

Abcb6 

AbcclO 

Abcc3 

Abcf3 

Abcg2 


Abhdl6a 

Abil 

Abll 

Abliml 

Abr 

Acaca 

Acadll 

Acatl 

AcotS 

Acp6 

Acsl4 

Actb 

ActrS 

Acypl 

Adaml7 

Adamtsl 

AdamtslO 

Adamts6 

Adamtsl4 

Adatl 

Adc 

AdcyS 

Adcy7 

Adk 

Ado 

Adora2b 

Adpgk 

Adrbkl 

Adss 


AF357374 

AF357376 

Afaplll 

Affl 

Aff4 

Afmid 

Aftph 

Aga 

AgbIS 

Agpatl 

Agt 

Ahcyl2 

Ahsa2 

AI118078 

AI450353 

AI462493 

AI597468 

Aicda 

Aifml 

AK006189 

AK006245 

AK007083 

AK013187 

AK015545 

AK016837 

AK018891 

AK018967 

AK019250 

AK019365 


AK031165 

AK035829 

AK037159 

AK038627 

AK040752 

AK042136 

AK043789 

AK043804 

AK043846 

AK043958 

AK044354 

AK044623 

AK045700 

AK047520 

AK048941 

AK051804 

AK053136 

AK053772 

AK054042 

AK078466 

AK079699 

AK079730 

AK079777 

AK085438 

AK087382 

AK089118 

AK131831 

AK132720 

AK134933 


AK140370 

AK141659 

AK142999 

AK154275 

AK155149 

AK155592 

AK156477 

AK158379 

AK162774 

AK163160 

AK164218 

AK165329 

AK166079 

AK167137 

AK172386 

AK188991 

AK204212 

AK212603 

AK212710 

Ak3 

Akapl 

Akapl3 

Akirinl 

Akirin2 

Akt2 

Aldoa 

Alg5 

AlkbhS 

Aloxe3 


Ambral 

Ammecrl 

Anapcl 

AnapcS 

Ang 

Ank3 

Ankrd24 

Ankrd40 

Ankrd52 

Ankrd54 

Antxrl 

Anxa2 

Anxa7 

Aplgl 

Apls2 

Ap3ml 

Ap4sl 

Aphlc 

Appll 

Arf3 

Arfgapl 

Arhgapl2 

ArhgaplS 

Arhgap21 

Arhgap26 

Arhgap4 

ArhgeflOl 

Arhgef25 

Arhgef? 


Atl3 

AtplBal 

Atplb2 

Atp2bl 

AtpSal 


Atp5b 

Atp6vOal 

Atp6vla 

Atp6vlb2 

Atp6vlel 

Atpaf2 

Atpifl 

Atraid 

Atrnll 

Atxn2 

Atxn3 

Atxn7ll 

Atxn7l3 

Avl9 

AW209491 

AW495222 

B230217C12Rik 

B230325K18Rik 

B330016D10Rik 

B3gnt2 

B930003M22Rik 

B930082K07Rik 

Bach2 

Banfl 

Baz2b 

BB019430 

Bbc3 

BC004004 

BC024582 


BC025920 

BC030867 

BC055111 

BC055823 

BC065397 

BC099561 

BC126883 

Bcarl 

Bcas2 

Bcas3 

Bcllla 

Bcl2ll 

Bcl7a 

Bdpl 

Bend3 

Bgn 

Birc6 

BIcap 

BImh 

Blocls4 

Bola2 

Bptf 

Bratl 

Brcal 

Brd3 

Brfl 

Bri3 

Brwd3 

Bscl2 


Btafl 

Btd 

Btf3l4 

C030013C21Rik 

C030034L19Rik 

C030037D09Rik 

Clra 

Clrb 

Clrl 

C2 

C2cd2l 

C2cd5 

C330006A16Rik 

C330007P06Rik 

Cabinl 

Cacnalf 

Cad 

Cagel 

Calm2 

Camk2d 

CapnS 

Capn9 

Capnsl 

Carl2 

Carhspl 

Cask 

Caskin2 

CaspS 

Cat 


CatsperZ 

Catspergl 

Cbfb 

Cbrl 

Cbx4 

Cbx7 

Ccdclie 

Ccdcl36 

Ccdcl4 

Ccdcl46 

Ccdcl76 

Ccdc22 

Ccdc65 

Ccdc77 

Ccdc84 

Ccdc92 

Ccdc93 

Ccnd2 

Ccnh 

Ccnk 

Ccser2 

Cct6a 

Cdl64 

Cdl64l2 

Cd247 

Cd276 

Cd2ap 

Cd300lh 

Cdcl4b 


Cdc25a 

Cdc26 

Cdc6 

Cdc73 

CdcaS 

Cdhl2 

Cdhl6 

Cdkllb 

Cdkl7 

Cdk2apl 

Cdk5rap3 

Cdkall 

Cdkl3 

CdkIS 

Cdknlb 

Cdkn2c 

Cdv3 

Celfl 

Cepl28 

Cfdpl 

Cfi 

Cflar 

Chd2 

Chd3 

Chekl 

Chfr 

Chracl 

Chrnal 

Chrnb4 


Chstll 

Chst3 

Chsyl 

Cic 

Cirbp 

Cisd3 

Cizl 

CkapS 

Claspl 

deal 

CIcfl 

Clcn2 

CIcnS 

Cldnl4 

Clic4 

Clintl 

Clipl 

Clk3 

CIptml 

CIspn 

CIta 

Cite 

Cmah 

Cmip 

Cnn2 

Cnot3 

Cnot6 

Cnot6l 

Cnppdl 


CnpyS 

Cnpy4 

Cog4 

Col8a2 

Commdl 

Commd? 

CommdS 

Comt 

Cops7a 

Cops8 

Copz2 

CoqlOa 

Coq5 

Cox7b 

CpebS 

Cpm 

Cptlc 

Cradd 

Creb3 

Crebbp 

Crebl2 

Creldl 

Crem 

Crk 

CrIfS 

Crnkll 

Croce 

Cryll 

Cryzll 


Csell 

CsgalnactZ 

Csnklal 

CsnklgS 

Csppl 

Csrnpl 

CtnnaB 

Ctnndl 

Ctsa 

Ctse 

CtsI 

Cttn 

Cuedcl 

Cuxl 

Cwcl5 

Cyb5 

CybSdl 

CybSrl 

Cyp2t4 

Cythl 

D030025P21Rik 

D030028A08Rik 

DMAbble 

D17Wsul04e 

D630003M21Rik 

D6Wsul63e 

D830025C05Rik 

Dapk2 

Dbi 


Dbp 

Dbrl 

Dcafl7 

Dcakd 

Dck 

Ddb2 

Ddit4 

Ddx20 

Debl 

Def8 

Denndlb 

Des 

Desil 

Dgatl 

Dgcrl4 

Dgka 

Dhodh 

DhrsB 

Dhx38 

Dhx57 

Dhx58 

Dhx8 

DIat 

Dleu2 

Dmtfl 

Dmtn 

Dnajbl4 

Dnajc25 

DnajcS 


DnajcS 

Dnml 

Dnmtl 

Dnttip2 

Dopey2 

Dpagtl 

Dpcd 

Dpepl 

Dpep2 

Dph3 

DpySO 

Dpysl2 

DpysIS 

DQ544183 

DQ548101 

DQ552992 

DQ553098 

DQ565977 

DQ566768 

DQ692659 

DQ694768 

DQ695092 

DQ695356 

DQ703023 

DQ704415 

DQ704975 

DQ705542 

DQ712916 

DQ713407 


DQ716966 

DQ725849 

Drl 

Drgl 

Drg2 

Drp2 

Dscr3 

Dtd2 

Dtnb 

Dtx3 

Dusll 

Dus2l 

Dus3l 

Duspl3 

Duspl6 

Dynlll 

Dynll2 

Dyrkla 

Dzip3 

E030042020Rik 

E130102H24Rik 

E130304l02Rik 

E2f2 

E530001K10Rik 

Eapp 

Ebfl 

Ecel 

Echl 

Ecsit 


Edeml 

Eed 

Eeflb2 

Eef2k 

Efcab2 

Ehbpl 

Ehmt2 

Eidl 

Eiflad 

Eif3c 

EifSd 

Eif3k 

Eif3l 

Eif4e 

Eif4ebp2 

Eif4g3 

Eif4h 

EifSa 

EifSb 

Elk3 

Elk4 

Ell 

Elmsanl 

Empl 

EntpdS 

Epb4.1 

Epb4.1ll 

Epcl 

Ephb4 


Epo 

Erbb2 

Ercc6l 

Ergicl 

Erlin2 

Ernl 

Esyt2 

Ets2 

Etv5 

Evi5 

Ewsrl 

Exd2 

ExoscS 

Ext  1 2 

Eya3 

F420014N23Rik 

Famllla 

Famll7a 

Faml29b 

Faml31a 

Faml34a 

Faml36a 

FamlSSa 

FamlSSa 

Faml9a2 

Fam214b 

Fam216a 

Fam35a 

Fam49b 


Fam57b 

Fam63b 

Fam73b 

Fam98a 

Fanca 

Fancc 

Fance 

Fancg 

Fars2 

Fbxll6 

Fbxl20 

Fbxo30 

Fbxo36 

Fbxo42 

Fbxo47 

Fbxo9 

Fchsd2 

Fdps 

Fes 

Fgd2 

Fgfrlop2 

Fhl3 

Fhl4 

Fig4 

Figf 

Fignll 

Fisl 

FkbplO 

Fkbpl4 


Fkbpla 

Fkbp5 

Fkbp7 

Fkbp8 

FIna 

Fmnl2 

Fnbpl 

FndcSa 

Fndc7 

Fnipl 

Fosl2 

Foxhl 

FoxjS 

Foxn2 

FoxnS 

Foxol 

Foxpl 

Fratl 

Frat2 

Frgl 

Fnnd4a 

Frs2 

Fry 

Fstll 

Fus 

Fut8 

Fxr2 

Fyttdl 

Fzd7 


G730013B05Rik 

Gabarapll 

Gabarapl2 

Gabpbl 

Gabpb2 

Gadd45b 

Gadd45g 

Gadd45gipl 

Gak 

Gapdh 

Gapvdl 

GaremI 

Gast 

Gata3 

Gatad2a 

Gate 

Gbas 

Gbp3 

Gcc2 

Gcic 

Gdf9 

Gdi2 

Gemin2 

Gemin6 

Gfpt2 

Gga2 

Ggnbp2 

Ggt5 

GimapS 


Ginsl 

Git2 

GIgl 

Gliprl 

Glol 

Gludl 

Gml0433 

Gm  10463 

Gml0610 

Gm  10642 

Gml0653 

Gml0655 

Gml0657 

Gml0658 

Gml0762 

Gml0837 

Gmlll84 

Gmll206 

Gmll292 

Gmll335 

Gmll336 

Gmll437 

Gmll453 

Gmll464 

Gm  11474 

Gmll491 

Gmll521 

Gmll602 

Gmll612 


Gmll619 

Gmll627 

Gmll630 

Gmll680 

Gmll696 

Gmll715 

Gmll827 

Gml2035 

Gml2054 

Gml2057 

Gml2060 

Gml2063 

Gml2245 

Gml2257 

Gml2279 

Gml2308 

Gml2309 

Gml2314 

Gml2358 

Gml2396 

Gml2795 

Gml2951 

Gml2974 

Gml2981 

Gml3054 

Gml3182 

Gml3201 

Gml3256 

Gml3297 


Gml3334 
Gml3363 
Gml3375 
Gml3398 
Gm  13447 
Gml3548 
Gml3559 
Gml3564 
Gml3626 
Gml3630 
Gml3657 
Gml3705 
Gml3770 
Gml3830 
Gml3836 
Gml3855 
Gml3936 
Gm  14005 
Gml4167 
Gml4216 
Gm  14455 
Gml4634 
Gml5411 
Gml5420 
Gm  15454 
Gml5688 
Gm  15747 
Gml5760 
Gml5787 


Gml5831 
Gml5860 
Gml5892 
Gml5903 
Gml5927 
Gml5962 
Gml6023 
Gml6185 
Gml6196 
Gml6197 
Gml6230 
Gml6274 
Gml6540 
Gml6557 
Gml6580 
Gm  16740 
Gml6880 
Gm  17077 
Gm  17098 
Gml7112 
Gml7138 
Gml7157 
Gml7300 
Gml7617 
Gml7661 
Gm  17705 
Gm  19705 
Gm20257 
Gm  20748 


Gm4221 

Gm4673 

Gm4978 

Gm5069 

Gm5134 

Gm5258 

Gm5428 

Gm5432 

Gm5464 

Gm5512 

Gm608 

Gm6225 

Gm6297 

Gm6444 

Gm6471 

Gm6525 

Gm7598 

Gm9812 

Gm9850 

Gm9900 

Gm9959 

Gm9985 

Gnal3 

Gnb2 

Gnb2ll 

Gnl3 

Golga3 

Got2 

Gpbarl 


Gpn3 

Gprl9 

Gpr35 

Gpr82 

Gpr85 

Gramdla 

Grb2 

GrcclO 

Grik4 

Gsel 

Gskip 

Gstt3 

Gtf2h2 

Gtf2i 

Gtf2ird2 

Gyk 

Gypc 

Gysl 

Gzmm 

HlfO 

H2afz 

H2-D1 

H2-DMbl 

H2-DMb2 

H2-L 

Hacll 

Hbpl 

Hdac7 

Hdgfrp2 


HeatrSa 

Helb 

Helz 

Herc4 

Hesl 

Hes7 

Hexim2 

Hid 

Hic2 

Hifla 

Hint2 

Hiplr 

Hira 

HiripS 

Histlh2ac 

Histlh2af 

Histlh2bb 

Histlh2bc 

HistlhSc 

Hivepl 

Hies 

HIx 

Hmbs 

Hmg20b 

Hmox2 

Hnrnpf 

Hnrnphl 

HnrnphS 

Hnrnpk 


Hnrnpl 

Hnrnpu 

Hnrnpull 

Hnrnpul2 

Homerl 

Hoxb6 

HoxbS 

Hoxc5 

HoxdlO 

Hoxd3 

Hps3 

Hsdl7bl2 

Hsp90abl 

Hspal3 

Hspa4 

Hspb9 

l830077J02Rik 

leal 

Icaml 

Id3 

Ifi35 

IftSO 

Igflr 

Igfbp6 

Ik 

lllrap 

Il3ra 

Ilf2 

Ilf3 


Ilk 

Immt 

Impdhl 

Ingl 

Ing3 

Ing4 

Inip 

InoSO 

InoSOd 

InoSOe 

InppSb 

IntsS 

Ipoll 

Iqce 

Iqcg 

Irf2bp2 

Irf2bpl 

Irs3 

Itfg2 

Itgb2 

ItgbS 

Itm2b 

Itpkb 

lzumo4 

Jarid2 

Kairn 

Kansll 

Kars 

Kat2a 


Kat2b 
Kbtbd? 
Kcnh3 
Kctdl9 
Kctd20 
Kdm3a 
Kdm4d 
KdmSa 
Kifll 
Kif23 
Kif24 
KifSb 
Kiel 
KlhdclO 
Klhlll 
KlhllS 
Kmt2ri 


Kntcl 


Let  I 

Ldbl 

Leol 

Leprel4 

Leprotll 

Let  m2 

Lhb 

Lias 

Lipe 

Lmf2 

Lmna 

Lmo2 

Lmo4 

Lmtk2 

Lnx2 

LOC100504608 

LOC100504703 

Loxl3 

LparS 

Lphnl 

Lpin2 

Lpxn 

Lrchl 

Lrgl 

Lng2 

Lrp2bp 

Lrrcl 

Lrrcl6a 

Lrrcl6b 


Lrrc46 

Lrrc49 

Lrrc58 

Lrrkl 

Lrsaml 

Lsmdl 

Luc7l 

Luc7l2 

Luc7l3 

Luzpl 

Ly6g6f 

Lynn4 

Lynn7 

LysmdB 

Lyzl6 

Mad2ll 

Mad2llbp 

Madd 

Maea 

Magtl 

Malatl 

Maml2 

Manea 

Map2k6 

Map2k7 

Map3kl2 

Map3k3 

Map4kl 

Map4k2 


Map4k3 

Mapk6 

Mapkbpl 

Mapt 

Marchs 

Marchs 

March? 

Marchs 

Mark? 

Marveldl 

Mast4 

Mbnl2 

Mbnl3 

Mbtps2 

Mcccl 

Mcm3 

Mcm3ap 

Mcm9 

Mctsl 

Mdpl 

Medl3 

Medl4 

MedlS 

Med6 

Melk 

Memol 

Mettll 

Mettll4 

Mettll? 


MettlS 

Mga 

Mical3 

Midn 

Mif4gd 

Mirl25a 

Mirl32 

Mirl52 

MirlSb 

Mirl7 

Mirl7hg 

MirlS 

Mirl931 

Mirl956 

Mirl99b 

Mirl9a 

Mirl9b-1 

Mir20a 

Mir212 

Mir3058 

Mir3109 

Mir5122 

Mir5135 

Mir615 

Mir670 

Mir677 

Mir702 

Mir92-1 

Mir99b 


Mirlet7e 

Mitdl 

Mkks 

Mlec 

Mlf2 

MlltlO 

MIxip 

Mmp25 

Mmsl9 

Mobla 

MobBa 

MorcB 

Morf4ll 

Morf4l2 

Mornl 

Morn2 

MornB 

Mpnd 

Mpp6 

Mpvl7 

Mrc2 
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